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PURPOSE 

One  example  of  a  SONAR  equation,  for  the  passive  detection  of  a  target 
which  is  radiating  sound  into  the  ocean,  is 


EL  -  TL  >,  NSL  +  10  tog  w  -  DI  +  DT 
where 

EL  =  Source  Level  of  the  target  being  detected  passively 

TL  =  Transmission  Loss  as  the  signal  propagates  to  the  detector 

NSL  =  Noise  Spectrum  Level  of  the  ambient  noise  in  the  ocean 

w   =  Frequency  Bandwidth  of  the  detecting  system 

DI  =  Directivity  Index  of  the  detecting  system 

DT  =  Detection  Threshold  of  the  detecting  system  for  a  specified 
probability  of  detection  and  a  specified  false-alarm  rate 

The  purpose  of  this  course  is  to  provide  the  physical  and  conceptual 
background  necessary  to  understand  the  meaning  of  each  of  these  terms 
and  how  they  are   combined  to  form  the  SONAR  equations.  At  the 
end  of  this  course,  the  student  should  be  able  to  apply  the  appropriate 
SONAR  equation  to  a  given  problem  and  estimate  such  things  as  maximum 
detection  range,  optimum  receiver  depth,  necessary  bandwidth,  etc. 
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where  *\r(^)    is  the  instantaneous  voltage  at  output 
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BEARING-RANGE  RECORDER  (active  sonar) 

brightness  related  to  the  strength  of 
the  echo  returned  from  each  range  and 
bearing 


LOFAR-GRAM  (passive  sonar) 

paper  is  electrically  burnt,  with 
darkness  related  to  the  strength  of 
the  signal  received  at  each  frequency 
and  time 


MONOFREQUENCY  SIGNALS 

For  ease  of  discussion,  we  will  at  first  restrict  ourselves  to  consideration 
of  tonal s  (mono frequency  signals).  Later  on,  we  will  see  that  more  compli- 
cated signals  can  be  broken  down  by  FFT,  Fourier  Analysis,  or  other  filtering 
techniques  into  a  collection  of  monofrequency  signals. 

The  time  history  of  the  acoustic  pressure  for  a  monofrequency  signal  can  be 
represented  as 
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P  =  The  effective  amplitude  of  the  acoustic  pressure.  This  is  not  the 
same  as  the  peak  amplitude  A,  but  is  more  convenient.  We  will  never 
use  the  peak  amplitude,  but  it  turns  out  that  P  =  0.707A  for  the  above 
signal. 

T  =  The  period  of  the  monofrequency  signal.  Tne  period  is  the  time 

interval  required  for  the  acoustical  signal  to  go  through  one  complete 
cycle,  returning  to  the  configuration  it  had  at  the  beginning  of  the 
interval . 


PLANE  WAVES 
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If       x  =  constant 


(For  simplicity,  take  x  =  0.) 
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If   £  =  constant 


(For  simplicity/  take  t  =  0.) 


^    -A  **(£*) 


As  time  increases,  the  wave  moves  towards  larger  x, 


A  K4VS  II  SP1CS  AND  TIM 


Speed  of  sound  c  depends  on  Iwptritigt, 
depth,  and  salinity 

Frequency      f  •  1A 
Wavelength  •     X 
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INTENSITY 


Intensity  =  power  per  unit  area 

where,  "area"  is  the  area  of  a  "window"  perpendicular 
to  the  direction  of  travel  of  the  wave . 

In  terms  of  measureable  quantities , 


where   PC   =  specific  acoustic  impedance  of  the  fluid 


J 


P<r      =  1.5  x  106  Rayl 
PC   =  415  Rayl 


for  water 
for  air 


1  Rayl  =  1  kg/m  sec 


SIZES  OF  THINGS  IN  WATER 
Speed  of  sound,  frequency,  and  wavelength: 

The  speed  of  sound  in  sea  water  is  approximated  by 

3 
c  ~  1.5  x  10  m/sec  _+  a  few  per  cent. 

It  must  be  remembered  that  because  of  depth,  temperature,  salinity,  and 

contamination  there  may  be  variations  up  to  a  few  per  cent  in  the  exact 

value  of  c. 


Approximate  sizes 

of  wavelengths  of  mono frequency 

sounds  of  different 

frequencies  are: 

Frequency 

Wavelength 

10  Hz  =  10  cps 

150  m 

50 

30 

100 

15 

500 

3 

1   kHz  =  1000  cps 

1.5 

5 

0.3 

10 

0.15 

50 

0.03  =  3  cm 

1  m  =  1.09  yd 

Pressure  amplitudes: 

/r   =1  atmosphere  =14.7  lb/ in  at  the  ocean  surface 


6 

=  10  lib    (yb  =  microbar) 


=  10  yPa   (yPa  =  microPascal) 


(0  increases  by  about  1  atm  for  every  33  ft  of  depth  in  water. 

Representative  acoustic  pressures  generated  near  a  sonar  transmitter: 

P  ~  104  to  106  yb  or  0.01  to  1  atm. 
If  P  ~  (j  ,  then  cavitation  is  likely  and  performance  may  be  degraded. 
Representative  weak  acoustic  pressures  detectable  by  a  sonar  receiver: 
P  ~>    10"2  yb  or  10"8  atm. 
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SOUND  PRESSURE  LEVEL 


SPL     =     20  log(P/Pref)   re  Pref 


0.0002  yb 

(old  noise  measurements) 

rex 

=  lyb 

(old  conventional  standard) 

lyPa 

(new  Navy  standard) 

Examples 

If  P  =  1  atm  =  106yb,  then 
r 
194  dB  re  0.0002yb 

SPL  =  <  120  dB  re  lyb 

1  220  dB  re  lyPa 

If  P  =  10*4  atm  =  102yb,   then 

114  dB  re  0.0002yb 
5PL  ={    40  dB  re  lyb 
140  dB  re  lyPa 

If  P  =  10*" *  atm  3  10"2yb,   then 

34  dB  re  0.0002yb 

SPL  =  (  -40  dB  re  lyb 

60  dB  re  lyPa 
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CONVERSION  BETWEEN  SPL's  WITH  DIFFERENT  REFERENCE  PRESSURES 


For  a  given  pressure  amplitude 


SPL  re  lyPa     =     SPL  re  Ipb     +     100  dB 
SPL  re  0.0002ub     =     SPL  re  lyb     +     74  dB 


Notice 

0  dB  does  not  mean  the  absence  of  sound  .  but  msans  P  =  P 

ref 

ref 
then 

SPL     =     20  log  (P/Pref)      =     20     log(Pref/Pref)     =      20  log  1     =     0  dB  re  P^f 
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TRANSMISSION  LOSS 


SPL(ri) 


SPL(r2) 


SPL(ri) 


FC*0 


20  log     F 

ref 

P(r9) 

20  log    p-i- 
ref 


P(rO 

SPLC^)         ,     20     lOg       pgiy 
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SAMPLE  CALCULATION  OF  TRANSMISSION  LOSS 


If  SPL  at  1  meter  =  120  dB  re  lyb 
and  SPL  at  range  r  s  40  dB  re  lyb 
then     TL     =120     -     40  dB     =     80  dB 


TL     does  not  depend  on  the  reference  pressure.     Let's  convert  the  above  SPL's 
to  reference  lyPa  and  recalculate  the  TL 

Now  SPL  at  1  meter  =  220  dB  re  luPa 
and  SPL  at  range  r  =  140  dB  re  lyPa 
so  that  TL     =     220     -     140     =     80  dB 
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TRANSMISSION  LOSS  FROM  SPREADING 

The  total  transmission  loss  TL  can  be  considered  to  arise  from  two  different 

causes : 

1)  As  the  rays  of  sound  propagate  out  from  the  source  and  travel  through 

the  water,  they  will  bunch  together  or  spread  out  depending  on  the  properties 

of  the  speed  of  sound  profile  for  the  water.  The  transmission  loss  which  arises 

from  this  effect  is  termed 

TL   =  transmission  loss  from  spreading 
g 

Some  simple  examples  of  spreading  for  rays  traveling  straight  lines  are 


Spherical  spreading 

TL  =  20  log  r 
g 


\ 


^3 


Cylindrical  spreading 
TL  =  10  log  r 


No  spreading 

TL  =  0  dB 
g 


In  the  real  ocean,  the  rays  never  travel  in  straight  lines,  so  that  these 
above  simple  examples  are  unrealistic.  We  will  see,  however,  that  there  are 
many  situations  in  which  the  sound  will  spread  out  spherically  or  cylindrically 
over  large  portions  of  its  path,  and  the  above  simple  equations  will  appear 
in  the  complete  transmission  loss  equation. 
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SPHERICAL  AND  CYLINDRICAL  WAVES 


SPHERICAL 
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TRANSMISSION  LOSS  FROM  DISSIPATION 
2)  There  are  a  variety  of  other  physical  effects  termed  dissipative 
effects,  which  reduce  the  SPL  of  the  sounds  in  addition  to  the  spreading 
of  the  rays.  Since  the  spreading  of  rays  is  described  by  TL  ,  we  will  call 
the  portion  of  the  transmission  loss  which  results  from  these  other  effects 
TL   =  Transmission  Loss  from  dissipation 

X/ 

There  are  a  number  of  these  which  are  important  in  underwater  sound  propagation: 

(a)  Absorption.  As  sound  travels  through  the  water,  there  is  a 
conversion  of  the  acoustical  energy  into  thermal  energy  because 
of  effects  which  are  very  similar  to  friction  in  moving  machinery 
and  hysteresis  losses  in  transformers. 

(b)  Transmission  into  the  bottom:  When  a  sound  beam  strikes  the  ocean 
bottom,  part  of  the  beam  is  reflected  back  upward  in  the  water 
and  another  part  is  sent  down  into  the  bottom.  The  beam  which 

is  reflected  back  upward  in  the  water  contains  less  energy 
than     the  beam  which  traveled  down  to  the  bottom. 

(c)  Scattering.  The  presence  of  rough  boundaries  at  the  ocean  surface 
and  bottom  can  lead  to  a  scattering  of  the  sound  out  of  the  reflected 
beam.  The  presence  of  bubbles  and  other  impurities  in  the  water 
itself  can  also  scatter  sound  energy  out  of  the  path  of  the  ray. 

(d)  Leakage.  In  certain  kinds  of  sound  propagation  in  the  ocean, 
like  the  surface  layer,  the  rays  of  sound  can  "leak"  out  the 
bottom  of  the  layer  each  time  they  bend  down  to  graze  it.  This 
"leakage"  results  in  a  loss  of  energy  from  the  layer. 
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OTHtR  LOSS    MECHANISMS 
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TRANSMISSION  LOSS 


In  general  the  total  transmission  loss  TL  can  be  composed  of  two  terms, 


TL  =  TL   +  TLn 


where 


TL   =  transmission  loss  from  geometrical  spreading 


TL?  =  transmission  loss  from  dissipation 
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SPEED  OF  SOUND,   TEMPERATURE,   AID  DEPTH 
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These  curves  show  how  the  speed  of  sound  is  related  to  depth  and 
temperature.  The  effect  of  salinity  is  not  included,  so  these 
curves  are  appropriate  only  for  water  of  constant  salinity. 
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SPEED  OF   SOUND   FORMULA 


C   «  \HH°\.  1 
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*         Temperature  Correction 
T  in  C 


i  2  /  C  -  3  ^  ^  Salinity  Correction 

' ~    *       '  S  in  parts  per  thousand 

■+ 

O    O  \  1    7  Depth  Correction 

z  in  meters 


Y 


Very  accurate  formulas  for  the  speed  of  sound  as  a  function 
of  temperature/  salinity,  and  depth  are  very  involved  and 
difficult  to  calculate.   The  above  expression,  although 
not  highly  accurate,  is  sufficient  for  many  calculations 
and  is  relatively  easy  to  use.   It  is  accurate  to  within 
about  0.1  meter  per  second. 

For  a  discussion  of  many  of  the  more  precise  formulas  and 
the  discrepancies  between  them,  the  interested  reader  is 
referred  to  V.  A.  Del  Grosso,  J.  Acoust.  Soc.  Am. 56, 
1084-1091  (1974)  . 
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TYPICAL   SOUND  5PF£D  PROFILE 
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HOURLY 


PROFILE  VARIATIONS 
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Speed  of  Sound 
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RAYS 
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RAYS  IN  CONSTANT  GRADIENTS 


c( 


J> 


yv 
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Snell's  Law: 
Along  any  ray 
c(z)/cos  © 
is  a  constant. 


R  is  the 

radius  of  curvature 

of  the  ray 


CW-  c.+  ft 


t 


is  the  gradient 


The  gradient  g  can  be  either  positive  or  negative. 
For  a  speed  of  sound  profile  composed  of  layers  -within  each  of 
which  the  gradient  has  constant  value,  then  the  segment  of  the 
ray  within  each  layer  is  the  arc  of  a  circle  whose  center  lies 
at  the  depth  z  -  -cQ/g  below  the  top  of  the  layer  (where  the 
speed  of  sound  is  cQ  at  the  top  of  the  layer  and  g  is  the  value 
of  the  fradient  for  that  layer).  The  radius  of  the  circle  is 


R  = 


<3.  CCr<L    & 

The  ray  always  bends  toward  the  lower  speed  of  sound 
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CONSTRUCTION  OF  RAY  PATHS 


Bottom 
layer 


Line  of  centers  for  top  layer 
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ACOUSTICAL  RECIPROCITY 


If  the  direction  the  sound  travels  over  the  propagation  paths  is 
reversed  by  exchanging  the  positions  of  the  source  and  receiver,  and  the 
propagation  conditions  remain  the  same,  then  the  transmission  loss 
between  source  and  receiver  will  remain  the  same. 


Active 


As  a  consequence,  the  transmission  loss  suffered  by  the  signal  in 
reaching  the  target  from  an  active  source  will  be  identical  with  the 
transmission  loss  experienced  by  a  signal  of  the  same  frequency  (and 
duration)  radiated  from  the  target  when  it  reaches  a  receiver. 
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RAYS  AND  A  SURFACE  DUCT 
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LOSS  FROM  THE  SURFACE  DUCT 


When  a  ray  reflects  from  the  ocean  surface ,  some  of  the  energy 
is  scattered  into  different  directions  and  can  propagate  out  of  the 
surface  layer.   This  energy  is  thus  lost  from  the  channel,  which 
increases  the  transmission  loss.   Additionally,  energy  can  diffuse 
out  of  the  bottom  of  the  layer.   These  effects  can  be  described  in 
terms  of  the  "loss  per  bounce".   The  number  of  bounces  is  the  range 
divided  by  the  skip  distance  r  . 

An  empirical  formula,  resulting  from  the  analysis  of  a  number 
of  experiments,  has  been  obtained  by  W.  F.  Baker,  J.  Acoust.  Soc. 
am.  57,  1198-1200  (1975)  — 


b  =  loss  per  bounce  in  dB 
SS  =  Sea  State 


This  equation  is  valid  within  a  restricted  range  of  parameters , 
given  approximately  by 

2  <  SS  <  5 

3  kHz  <  f   <  8  kHz 

80  ft  <    layer  depth  <  200  ft 

30  ft  <  source  depth  <   60  ft 
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TRANSMISSION  LOSS  IN  THE  MIXED  LAYER 


K 


^ —  transition  range  r. 


D   r 


skip  distance   r 

s 


The  depth  of  the  mixed  layer  is  designated  by  D. 

The  depth  of  the  source  in  the  layer  is  designated  by  z  . 

The  transition  range  r  is  the  distance  the  sound  beam  must  travel 
until  it  effectively    fills  the  channel.  It  is  calculated  from 

A  -  \,0'5*\0l\\   — — (a11  distances  are  in  meters) 

t        V  T>-Ji  - 

The  skip  distance  r  is  the  distance  required  for  the  rays  which  just 
graze  the  bottom   '  of  the  layer  to  complete  one  full  cycle  from 
the  layer  bottom  to  the  surface  and  back  to  the  bottom.  It  is 
calculated  from 

/I  -  2.H  ffO    J'j)  (a11  distances  are  in  meters) 

If  the  loss  per  bounce  is  designated  by  b  and  the  energy  loss  resulting 
from  absorption  is  designated  by  a,  then  the  transmission  loss  for  the 
mixed  layer  can  be  estimated  by  tie  following  formulas: 


TL  *  20  log  r  +  ar  within  transition  range 

TL  «  10  log  r   +  10  log  r  +  ar  +  br/r     beyond  transition  range 


For  ranges  less  than  r  the  sound  is  spreading  out  spherically  to  fill 
the  sound  channel.      For  ranges  greater  than  r  ,  however,  the  sound 
has  spread  far  enough  vertically  to  be  reflected  ^  from  the  surface  and 
bent  back  upward  from  the  bottom  of  the  layer.  Once  this  hap  ens,  then 
the  sound  is  trapped  and  will  spread  cylindrically. 
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SAMPLE  CALCULATION   OF  A  TRAESMISSIOE    LOFfT 


Assumptions: 


Depth  of  the  mixed  layer  =  i^O  n 
Depth  of  the  source      *  20  m 
Frequency  cf  the  sound   =  2  kHz 
Sea  state  ■  3 

1  e  then  determine 

Skip  distance       =  5300  m 

Transition  ran£e     ■  94O  m 

attenuation  constant  -  0.00015  dB/m 

Loss  per  bounce      »=  3.5  dB/bounce 

The  transmission  loss  formulas  are  thus 

(a)  For  ranges  less  than  the  transition  ran^e, 

TL  -  20  log  r  +   (1.5  x  10'u)r 

(b)  For  ranges  greater  than  the  transition  rar.ge, 

TL  *  10  log  r  +  (8.1  x  10-k)r  +  30   . 
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TRANSMISSION  LOSE  CURVE  FOR  SOURCE  AND  RECEIVER  IN  TIE  DUCT 


Parameters: 


Depth  of  the  mixed  layer  ■  1*0  m 
Depth  of  the  source      ■  20  m 
Frequency  of  the  sound   ■  2  kHz 
Sea  state  =  3 

Receiver  is  in  the  duct 


JOO 

7  ) 


7 


Lo 


1006 
i —  i      i '  i"  r>  — 


10,00  0 


znrit     (tr) 
> 


Simple  model 


to 


to  u    - 


Empirical  curve  representing 
the  average  of  vast  amounts  of  data,  from 
Navy  Underwater  Sound  Lab.  Tech.  Mem.  1110-14-55,  1955 


BUT:  These  results  can  be  radically  changed  by  effects  such  as 
Surface  Interference  (which  becomes  important  at  lower 
frequencies  and  lower  sea  states),  steep  gradients  below 
the  layer,  and  source  and/or  receiver  near  the  top  or 
bottom  of  the  duct.  Additionally,  if  the  frequency  of 
the  sound  is  too  low,  it  will  not  be  trapped  in  the  duct. 

As  a  very  rough  estimate,  the  frequency  must  be  greater  than  about 
(2  x  1()5)/d3/2,  where  D  is  in  meters  and  the  frequency  is  in  Hz, 
for  the  sound  to  be  well  trapped  in  the  layer  and  for  the  above 
model  to  apply. 
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SURFACE  INTERFERENCE 


m 


<&"■ 


f* 


r 


SOURCE 


<a^  ^    > 


,S 


o^l 


'4  I 


RlCtlVBf^ 


*<- 


r 


f?     ~      sovRea    tit  p  i*i 


d      ~  R£C£li/&C      pjUL'-T't 


a? /^ A//?z--     f  yoKaoN 7 At  ) 


34 


Acoustic  pressure  at  the  receiver 
P 


direct  ray- 


reflected  ray 


T'/*U 


pressure  amplitude  *  A 


Ad  >  Ar 

Ad  *Ar 
Ad  »  Ar 


because  (1)  reflected  rays  travels  a  longer  path 

(2)  scattering  from  rough  surface 
if  range  large  and  surface  smooth 
if  surface  very  rough 


phase  delay  -  ^   (due  to  difference  in  path  length  and  to  phase  change 

upon  reflection) 


A  -  oc 


A    -  l8(f 
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rteal  sound  speed  profiles  will  affect  TL  only  at  large 
distances,      (i-otice  the   shadow  zone  predicted  at  6"  km.) 
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The  range  at  vihich  the  received  signal  is  a  minimum  is  £iven  by 

zhJ 


-n  c 


J  vhere      ty    '   '  ?   **■  ,    J  , 


It  is  pos?sible  to  use  information  from  the  received  signal  to  determine  the 
depth,  speed,  and  distance  of  closest  approach  of  the  source.  For  a  source 
moving  at  constant  speed  v  along  a  strai£ht  line  at  a  fixed  depth  h,  the 
rar:{_e  r  is  given  by 


ra  r   j*  i  r>j-xt 


vhere  at  t  =  0,  r  ■  s  -  the  distance  of  closest  approach. 

Combining  these  equations  ^ives  the  time  at  which  the  signal  is  a  minimum  as  a 

function  of  frequency. 
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The  display  that  results  is  shown  below: 
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Our  equation  can  be  rewritten  in  standard  form 


-  / 


vhere 


Q     r 


and 


SA 


ir 


3  hd 

This  is  the  equation  of  a  hyperbola 


vhere  a  is  the  value  of  f  at  t  -  0  and  t/a  is  the  slope  of  the  asymptote, 

If  v  is  calculated  from  the  measured  Doppler  shifts  of  the  discrete  lines,  * 

the  depth  of  the  source  is  found  from 

I 


'  ■   *  (If)  i 

and  the  distance  of  closest  approach  from 


S   - 


CL 


Cf  course,   for  day-to-day  use,   nomograms   can  be  developed  to  simplify  the 
calculations. 


*      It  will  be   shown   later   that 

m  -  2  r  f 

where  &$     is  the  total  frequency  change  of  a  line. 
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PROPAGATION   PATHS  FOR  SHALLOW  SOURCE   Ai\D  RECEIVER 


c(z) 


Range 


(a)  Surfade  Duct.  The  isothermal  water  at  the  ocean  surface  forms 
the  mixed  layer  in  which  sound  can  be  trapped.  The  continual 
reflection  at  the  rough  surface,  refraction  at  the  bottom  of 
the  layer,  and  leakage  provide  relatively  large  energy  losses 
so  that  this  path  is  useful  for  relatively  short  ranges. 

(b)  Convergence  Zone.  Vhen  the  water  is  deep  enough,  the  sound  rays 
■which  travel  down  past  the  SOFAR  channel  axis  can  be  bent  back 
upward  before  they  reach  the  bottom.  These  rays,  when  they 
approach  the  surface  again,  tend  to  be  brought  together  so  that 
the  sound  energy  becomes  more  concentrated,  and  the  transmission 
loss  is  reduced  by  a  significent  amount.  Under  advantageous 
conditions,  these  rays  will  be  refocused  in  second,  third,  and 
higher  convergence  zones. 

(c)  Bottom  Bounce .  The  reflection  of  sound  off  the  bottom  and  back 
to  the  surface  provides  a  propagation  path  filling  in  the  ranges 
between  the  greatest  useful  range  of  the  mixed  layer  and  the 
innermost  edge  of  the  convergence  zone.  Unfortunately,  the 
unpredictable  and  usually  large  reflection  loss  that  occurs  at 
the  bottom  causes  the  transmission  loss  along  this  path  to  be 
rather  large. 


40 


/  .1 

fir* 


75   - 


2> 


rs 


CO 

w 
o 
►-} 

c 
o 

•H 
G; 
(0 


c 


O-i 


/dfl 


J/0 


f2a 


I  -  2  kHz  a  -  0.08  dBAm 

Vater  depth  «  15,000  ft  -  i;600  m 

Layer  depth  «  100  ft  -  30  m 

Source  and  receiver  depth  «  50  ft  ■  1$  m 

Bottom  loss  ■  5  dB/bounce 
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\   Spherical  spreading  alone 
\ 
\ 
\ 
\ 


\  Spherical  spreading 


plus  absorption 
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),ooo 


Range  (Ion) 


(a)  Surface  duct 

(b)  Convergence  zone 

(c)  Bottom  bounce 


Ref .     Urick,    Principles  of  Underwater  Sound  for  Engineers 
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CQNVERCEFCJS  ZONE  PARAMETERS 


\<-AC  -X 


> 


c  ■  the  maximum  speed  of  sound  above  the  SOFAR.  channel  axis 

A  c  »  the  difference  between  c  and  the  speed  of  sound  at  the  elbow 

Z  ■  the  vertical  distance  between  depths  for  which  c(z)  =  c 

r   »  the  ran^e  to  the  conver£ence  zone 

Snell's  Law  reveals  that  a  ray  which  is  horizontal  in  the  upper  ocean  where 
c(z)  ■  c  will  also  be  horizontal  at  a  distance  Z  deeper,  and  a  ray  which 
is  directed  downward  where  c(z)  *  c  must  penetrate  below  the  distance  Z. 
To  insure  that  enough  rays  are  returned  to  the  surface  to  focus  at  the 
convergence  zone,  the  total  depth  of  water  must  exceed  Z  or  Z  +  D  by  an 
amount  called  the  depth  excess ^  d 

If  there  is  no  surface  layer,  then  c  is  associated  with  the  surface 
and  Z  is  measured  from  the  surface.  If  there  is  a  surface  layer,  then  c 
is  associated  with  the  bottom  of  the  layer  and  Z  measured  from  the 
bottom  of  the  layer. 


r«-AC  -H 


Z+P- 


-777777-7777T  ~  ~ 
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COKVEBGERCE  ZONE  RAKGjg!  Al'D  TRANSMISSION  LOSS 


A  rough  estimate  of  the  range  at  which  a  convergence  zone  forms  is  given 


r  „  ++ 

CZ 


.8  Z  \Jc  /Ac       +   Layer  effect 


Layer  effect  ■  0  if  there  is  no  surface  duct 


r  if  a  surface  duct  is  present  and  both  source 
and  receiver  are  near  the  top  of  the  layer. 


The  accuracy  of  the  above  expression  depends  on  how  veil  tue  constant- 
gradient  approximations  fit  the  true  c(z)  profile.  As  might  be  expected, 
when  the  elbow  of  the  curve  is  very  sharp,  and  the  gradients  above  and 
below  the  elbow  are  relatively  constant,  then  the  prediction  will  be 
rather  good.  The  width  of  the  convergence  zone,  over  which  the  rays 
are  appreciably  focused,  is  normally  about  10  per  cent  of  the  range  to 
the  zone. 

Since  the  rays  begin  to  spread  out  spherically  as  they  leave  the 
source  and  are  attenuated  only  by  the  absorptive  losses  in  the  ocean, 
the  transmission  loss  for  convergence  zone  paths  can  be  written  in  the 
form 


TL  «  20  log  r  t  ar  -  G 


where  C  is  the  convergence  zone  gain.  Typical  values  of  G  range  between 
about  $   and  l£  db  when  a  convergence  zone  is  found. 
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Typical  Convergence  Zone  Ranges 


North  Pacific 


Surface  Temperature-: 


50°  F 

55 

60 

65 

70 

75 

80 


Minimum  depth  of  CZ  operation* 


1270  ftm  «  23? 0  m 

1610  2 940 

1900  3R80 

2150  3930 

2R00  i;uOO 

2600  U760 

2800  5130 


Range 

47  kyd  ■ 

1 
to  first  CZ 

=  43  km 

52 

47 

56 

51 

60 

^ 

64 

57 

66 

60 

69 

63 

Area 


1!.  Pacific 

! 

N.  Atlantic   ( 

I 
Norwegian  Seal 

^Mediterranean 


Surface  lemperature-: 
50°  F 

! 

50 
50 

67 


•linimum  depth  for  CZ    Range  to  first  CZ 


operation* 
1270  ftm 
1050 
1680 

800 


47  kyd 
46 
S3 
33 


J, 


*  Assumes  surface  duct  absent 

+  Allows  for  a  200  ftm  depth  excess 

Ref .  Convergence  Zone  Range  Slide  Pule,  Naval  Indervater  Systems  °enter 
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BOTTOM  POUNCE 


For  ari£les  of  depression  beyond  about  15>   ,   the  radius  of  curvature  of 
a  ray  is  so  lar^e   that  it  is  a  ^ood  approximation  to  assume  straight-line 
propagation  of  the  rays  over  the  distances  of  interest. 


A  ray  bouncin£  off  the  bottom  -will  return  tc  the  surface  at  a  ran^e 

r  «=  2   H  tan  6 

Since  the  distance  travelled  by  the  ray  is  r/(cos  ©),  and  since  the 
rays  spread  out  spherically  and  are  attenuated  only  by  the  absorptive 
losses  of  the  ocean  and  the  reflection  at  the  bottom,  an  approximate 
transmission  loss  can  be  written  as 


TL  =  20  lp£  — JL-  +  —^   +  BL 

cos  0     cos  0 


vhere  BL  is  the  bottom  loss  in  dB/tounce. 
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BOTTOM  LOSS 


Depen  is  on  frequency  and  £  razing    an^ le 
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Depends  on  bottom  composition 

Example:  At  2\\   kHz  and  10°  grazing  an£le 


Mud 

16  di 

./bounce 

Mud -sand 

10 

Sand -mud 

6 

Sand 

h 

Stony 

k 

Depends  on  surface  roughness 

Bottom  loss  [.oes  up  as  ratio  of  roufhness-to-viavelength  increases 


Kef.  Urick,  Princioles  of  Underwater  Sound  for  Engineers 
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RSUABLE  PROPAGATION  PATH  FOR  DEEP  SOURCE  OR  RECEIVER 


c(z) 


Range 


t 


Reliable  Acoustical  Path  (RAP). 

If  the  target  is  shallow  and  the  sensor  deep,  or  vice  versa,  there 
exists  the  possibility  of  a  guarenteed  acoustical  path  between  the  two. 
This  path  is  useful  to  ranges  beyond  those  attainable  by  transmission  in 
the  surface  duct.  As  the  deeper  of  the  two  (sensor  or  target)  is  lowered 
further  and  further  below  the  layer,  the  shadow  zone  indicated  in  the 
diagram  moves  out  to  larger  ranges.  If  there  is  no  mixed  layer,  then  the 
leading  edge  of  the  shadow  zone  is  formed  by  the  downward-bending  ray 
which  just  grazes  the  surface. 

Because  the  angles  for  RAP  are  usually  fairly  large,  the  rays  can 
be  considered  to  travel  approximately  in  straight  lines  between  the 
target  and  sensor.  For  this  kind  of  propagation,  the  rays  spread  out 
spherically  and  the  only  losses  are  those  from  the  absorption  of  sound 
in  sea  water.  An  approximate  transmission  loss  formula  is  then 


TL 


20  log  _r    ♦  JSL_ 
cos  ©    cos  © 


where  0  is  the  angle  of  elevation  of  the  rays  between  target  and  sensor. 

If  the  negative  gradient  (below  the  layer,  if  one  exists)  has 
magnitude  g,  then  the  range  to  the  shadow  zone  at  the  surface  is  estimated 

by  . , 

1.U  Vc0(z-D)/g   +  h  rs 

where  (z  -  D)  is  the  depth  below  the  bottom  of  the  layer.  If  there  is  no 
layer,  then  D  «  0  and  there  is  no  skip  distance,  so  r  ■  0, 
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THE  SOFAR  CHANNEL 


c(z) 


Range 


Just  as  for  the  surface  duct,  the  sound  spreads  out  more  or  less 
spherically  until  the  rays  are  bent  back  into  the  channel,  after  which 
the  sound  spreads  cylindrically.  Thus,  beyond  some  transition  range  r.  ' 
we  should  have  a  spreading  term  -which  goes  as  10  log  r.  Those  rays 
whose  angles  of  inclination  with  respect  to  the  SOFAR  channel  axis  are 
large  enough  that  they  reflect  from  the  ocean  surface  or  bottom  will 
suffer  reflection  losses  and  scattering  losses  and  eventually  be  so 
weak  compared  to  the  trapped  rays  that  they  can  be  neglected.  Thus, 
for  propagation  in  the  channel,  there  are  no  losses  except  absorption, 
so  that  the  transmission  loss  has  the  simple  form 


TL  =  10  log  r  +  10  log  rt« 


ar 


for  ranges  beyond  r.  ' 


The  transition  range  depends  on  how  deep  the  sound  source  is  in  the  channel 
and  the  vertical  extent  of  the  channel.  For  sources  placed  reasonably 
down  in  the  channel,  the  transition  range  is  less  than  the  convergence 
zone  range  r  .  Notice  that  the  rays  cycle  up  and  down  in  the  channel, 
accomplishing  one  cycle  over  an  interval  very  similar  to  the  convergence 
zone  range . 


1  aa*  & 
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TIMS  STRiiTCmG  III  TIE  S6FAR  CilANhEL 


L 


^Ni 


time 


An  interesting  fact  often  observed  in  the  SOFAR  channel  is  that  transient 
signals  are  "stretched"  in  time  as  they  propagate  to  large  ranges.  This  arises 
because  of  the  variation  in  the  speed  of  sound  over  the  vertical  extent  of  the 
channel.  For  the  SGFAfr  channels  commonly  encountered  in  the  northern  hemisphere, 
the  speed  of  sound  increases  sufficiently  rapidly  with  distance  above  or  below 
the  channel  axis  that  sound  traveling  otev   those  rays  which  swing  over  large 
vertical  distances,  like  a  and  b  in  the  above  sketch,  reaches  the  receiver 
at  times  earlier  than  sound  traveling  over  the  more  direct  paths  c,  d,  and  e. 
Thus,  the  sound  that  travels  straight  down  the  axis  of  the  channel  arrives  at 
the  receiver  last.  This  leads  to  the  characteristic  blossoming  of  the  received 
signal  from  an  explosive  charge,  the  sound  first  being  very  weak,  and  then 
growing  in  volume  until  there  is  a  sudden  cutoff  of  sound  as  the  last  signals 
come  over  the  slowest  paths. 

The  exact  details  of  these  effects  depend  greatly  on  the  character  of  the 
speed  of  sound  profile:  If  the  sound  does  not  increase  sufficiently  rapidly 
above  and  below  the  axis,  for  example,  the  time  stretching  may  not  occur,  or 
may  even  be  reversed,  with  the  sound  traveling  straight  down  the  axis  arriving 
first.  Additionally,  the  exact  location  of  the  source  and  receiver  with  respect 
to  the  channel  axis  is  very  important  with  respect  to  the  envelope  of  the 
received  signal  and  the  amount  of  time  stretching  encountered.  For  the 
typical  profiles  in  the  northern  hemisphere,  the  increased  duration  At  of  the 
received  signal  over  that  sent  can  be  estimated  by 


At 
t 


1  A_c 

3  c 


if  source  and  receiver  are  on  the  channel  axis,  and  less  if  they  are  not. 
The  quantity  t  is  the  total  time  of  flight  between  source  and  receiver  and 
Ac  is  the  maximum  variation  in  the  speed  of  sound  between  the  upper  and 
lower  boundaries  of  the  channel  indicated  in  the  previous  page.  Typically, 
this  leads  to  about  °  seconds  of  stretching  for  each  1000  nautical  miles 
the  signal  has  travelled. 
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NORMAL  MODES 

Another  approach  to  the  propagation  of  sound  in  the  ocean  is 
Normal  Mode  Theory.  If  the  mathematical  difficulties  of  this 
approach  can  be  overcome,  it  yields  more  exact  values  for  the 
transmission  loss  than  can  be  obtained  using  Ray  Theory. 

Normal  Mode  Theory  yields  a  sum  of  individual  solutions 
(each  called  a  normal  mode)  which  must  be  combined  to  obtain 
a  complete  expression  for  the  sound  field. 


•f  r  fi^fi 


«*  t  *   <j>   + 


f»+' 


hi 

=  Z-f* 


Each  one  of  these  normal  modes  has  its  own  behavior, 


A 


f„-Vf^^  [^a-vi 


^_ y . 

This  part  looks  just  like  a  cylindrical  wave 
traveling  outward  with  its  amplitude  falling 
off  like  ^/\/f%  .  Notice  however  that  it  has 
its  own  special  speed  cn. 

This  part  is  an  amplitude  factor  which  changes  with 
depth.  Notice  that  each  normal  mode  has  its  own 
special  depth-dependent  amplitude  factor. 


These  special  factors  Zn  are  called  "normalized  eigenf unctions" 
They  are  found  by  solving  the  wave  equation  and  specifying  the 
behavior  of  the  pressure  at  the  ocean  surface  and  bottom.  The 


solution  for  each  Z  also  gives  a  required  value  for  c 


the  depth  dependence  of  Z  and  the  value  of  cn  depend 
frequency  of  the  sound  waves. 


n* 


Both 


on  the 


This  kind  of  a  solution  can  be  obtained  whenever  the  speed  of 
sound  profile  c(z)  is  shaped  so  that  the  sound  can  be  trapped 
in  a  channel.  The  trapping  of  sound  in  the  mixed  layer  is  one 
such  example,  and  the  propagation  of  sound  to  long  ranges  in 
the  SOFAR  channel  is  another. 
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NORMAL  MODES  FOR  A  CHANNEL 
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The  above  sketches  represent  the  normalized  eigenfunctions  of  the 
first  three  and  sixth  normal  modes  for  sound  trapped  in  the 
SOFAR  channel.  These  sketches  show  that  the  greatest  part  of 
the  pressure  field  is  found  in  the  channel,  but  the  long  "tails" 
which  extent  to  depths  far  below  the  channel  axis  show  that 
there  is  sound  ensonification  at  depths  greater  than  those 
predicted  by  drawing  ray  paths. 
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CUTOFF  FREQUENCY  AND  LOSSES 
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Each  normal  mode  has  its  own  cutoff  frequency  f  .   If  the  sound 
has  a  frequency  below  that  of  cutoff  for  a  particular  mode,  then 
that  mode  cannot  be  excited  to  carry  sound  energy  in  the  channel. 

If  the  frequency  of  the  sound  is  slowly  increased  with  time,  the 
above  sketches  show  how  a  normal  mode  can  be  excited.  Just  above 
cutoff,  the  normal  mode  has  a  very  long  tail  which  goes  to  great 
depths  beyond  the  edges  of  the  channel.  As  the  frequency  increases 
further  above  cutoff,  the  tail  gets  shorter  and  shorter,  until 
at  frequencies  well  above  f  there  is  very  little  sound  existing 
outside  the  channel. 

This  means  that  for  a  specified  frequency,  sound  can  be  trapped  and 
sent  to  large  distances  only  in  those  normal  modes  whose  cutoff 
frequencies  are  less  than  the  frequency  of  the  sound.  Normal  modes 
with  higher  cutoff  frequencies  will  have  longer  tails  than  those 
with  lower  cutoff  frequencies,  and  those  whose  cutoff  frequencies 
exceed  the  frequency  of  the  sound  cannot  be  excited  at  all. 
The  larger  the  tail  extending  out  of  the  channel,  the  more  energy 
can  be  lost  from  the  sound  in  the  channel.  This  means  that 
higher  modes  (those  with  higher  cutoff  frequencies)  will  attenuate 
faster  and  be  less  significent  at  great  distances.  This  effect  is 
sometimes  called  mode  stripping. 

These  energy  losses  for  each  mode  can  be  handled  by  specifying  an 
absorption  coefficient  ^^  for  each  normal  mode.  Then,  we  have 

where  or_  increases  with  increasing  n. 
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TRANSMISSION  LOSS  (NORMAL  MODES) 
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RAYS,  NORMAL  MODES,  AND  REALITY 

Normal  mode  theory  is  a  complicated  mathematical  formalism  which  does  not 
lend  itself  to  simple  intuitive  interpretation.  However,  it  is  an  exact 
theory  which  avoids  certain  approximations  inherent  in  ray  theory.  It  is 
useful  in  discribing  situations  for  which  ray  theory  fails.  An  example 
of  one  such  situation  is  illustrated  below. 


Range 


The  source  is  in  the  layer  and  the  receiver  is  below  the  layer.  Ray 
theory  predicts  a  shadow  zone  at  large  enough  ranges.  While  there  can 
be  some  scattering  from  the  ocean  surface  into  the  shadow  zone,  this 
will  be  small  for  low  frequencies  and  low  sea  states.  The  surface 
interference  pattern  at  short  ranges  shows  that  the  sea  was  quite 
smooth  (the  observed  sea  state  was  about  1  or  2). 


Ray  Theory 

Normal  Mode  Theory 

Experimental  Data 


f  -  530  Hz 
D  -  300  ft 
-  200  Hz 


U00  ft 


too 


Range  in  kiloyards 


Notice  that  the  experimental  data  show  significent  ensonification  of  the 
shadow  zone.  This  cannot  be  predicted  by  the  ray  theory,  and  the  scattering 
of  sound  from  the  surface  is  not  strong  enough  to  account  for  it. 

Reference:  Pederson  and  Gordon,  Jour.  Acoust.  Soc.  Am.  37>  105>  (1965) 
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SUMMARY  OF  TRANSMISSION  LOSS  EXALTLI'CS 


For  both  source  and  receiver  in  the  surface   layer,   and  at  ranges 
greater  than  the  transition  range,  the  transmission  loss  can  be 

estimated  by 


TL     ■     10  log  r.      +     10  log  r     +     ar     +     br/r   . 

T>  S 


For  a  conver£ence  zone  detection,  we  can  vrite 
TL  ■=  20  log  r  +  ar  -  G  . 


If  the  propagation  path  is  bottom  bounce,  and  the  angle  of  depression 
of  the  sound  beam  in  that  mode  is  sufficiently  large,  then  the  simple 
geometrical  model  fives 

TL  -  20  log  -~g  +  — £~  +  BL  . 
L  cos  9     cos  9 


For  propagation  between  a  deep  point  and  a  shallow  point  over  the 
reliable  acoustical  path,  an  approximate  transmission  loss  formula  is 

TL  -  20  log  — L-,  + 


cos  9    cos  9 


Finally,  for  the  propagation  of  sound  trapped  in  the  S0FAR  channel, 
if  the  range  is  greater  than  the  transition  range  r. ',  the  transmission 
loss  can  be  estimated  by 


TL  -  10  log  r  +  10  log  rt'  +  ar 
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COMBINING  SIGNALS 

Now,  let  us  no  longer  restrict  ourselves  to  tonals  (monofrequency 
signals).   Because  we  really  don't  deal  with  only  a  single  frequency, 
we  have  to  find  out  how  signals  of  different  frequencies  combine. 
Consider  the  combination  of  a  large  number  of  monofrequency  signals t 


=  ^/v^Wv^ 


w  yyvyywwy>J 


The  combination  of  many  signals,  each  a  pure  tonal,  will  result  in 
a  total  signal  which  does  not  appear  to  have  definite  amplitude 
or  period.   However,  the  composite  signal  is  an  acoustical  wave, 
and  the  total  acoustical  energy  it  possesses  is  the  sum  of  the 
acoustical  energy   of  each  of  the  individual  tonals.   Thus, 
if  Em  is  the  total  energy  of  the  composite  signal,  we  have 


Em  =  El   +  E; 


M 

i'  I 


where  E^  is  the  energy  of  signal  1,  E2  that  of  signal  2,  and  so  on, 

For  the  complicated  signals  encountered  in  real  life,  which  have 
large  variations  from  moment  to  moment  in  their  amplitude  and 
shape,  the  frequencies  of  the  individual  tonals  which  combine  to 
form  them  are  extremely  closely  spaced.   Indeed,  it  is  often 
more  convenient  to  talk  about  the  energy  contained  in  the 
collection  of  tonals  which  lie  between  some  specific  frequency  and 
another  frequency  which  is  only  inf initisimally  different. 
If  aE  is  the  total  energy  of  all  the  tonals  lying  within  the 
incremental  bandwidth  aw,  then  we  can  define  the  quantity  e 
as  the  average  energy  density  over  the  incremental  bandwidth  aw  by 


£  -- 


&JUJ 


Notice  that  the  energy  density  has  the  units  of  energy/Hz.   If  we 
analyze  some  complicated  signal,  therefore,  we  can  describe  it 
in  terms  of  the  energy  density  it  contains  in  each  incremental 
interval  of  frequency  over  the  entire  range  of  frequencies  which 
comprise  it. 
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ENERGY  SPECTRUM  AND  PRESSURE  SPECTRUM  LEVEL 

An  example  of  the  analysis  of  some  complex  signal  into  the 
energy  densities  found  in  each  frequency  interval  is  suggested 
in  the  following  curve i 
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The  total  energy  of  the  signal  is  given  by 


ET   -   £,4>ur  +lxbMTx 


-f  *  •  • 


=  !*• 


&U4T. 


which  is  the  area  under  the  curve  of  energy  density  vs  frequency. 

If  all  of  the  ^w^  are  chosen  to  be  1  Hz  wide,  then  the  graph  gives 
the  energy  to  be  found  in  each  1  Hz  interval,  and  the  curve  is 
called  the  energy  spectrum  of  the  signal. 

It  is  more  conventional  in  acoustics  to  describe  the  signal  not 
in  terms  of  the  energy  spectrum,  but  rather  in  terms  of  the 
pressure  spectrum  level  PSL,  which  we  shall  define  here  as 


PSL  --  HO   A*£^* 


H.f 


JB 


re 


re 


re 


where  the  quantity  p  ,  the  effective  pressure  per  Hz,  generates 
the  associated  energy  found  in  the  1  Hz  bandwidth.   For  example, 
if  we  receive  a  certain  acoustical  signal,  pass  it  through  a 
filter  with  a  1  Hz  bandwidth,  and  obtain  the  effective  pressure 
amplitude  of  the  resultant  filtered  signal,  say  100  ^*.Pa,  then  the 
value  of  p   for  that  1  Hz  bandwidth  is  100  ^Pa/l  Hz  *  100  /*Pa/Hz. 
A  representative  curve  of  the  pressure  spectrum  level  vs  frequency 
is  presented  below: 


Tonal  or  "Spike" 


^0,. 


<XA 


— •  o 


Frequency 
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AVERAGED  ENERGY  SPECTRUM  AND  AVERAGED  PSL 


Very  often,  what  is  presented  as  the  energy  spectrum  or  the 
pressure  spectrum  level  is  not  that  which  would  be  measured  with 
a  1  Hz  filter,  but  is  rather  calculated  from  measurements  made 
with  a  filter  of  larger  bandwidth. 

The  overall  energy  to  be  found  in  some  arbitrary  bandwidth 
w  is  the  area  under  the  energy  spectrum  vs  frequency  curve 
over  the  appropriate  range  of  frequencies: 


+  jur->- 


Frequency 


This  gives  the  total  energy 


Ew  found  within  this  bandwidth,  and 
does  not  say  how  E   is  distributed  over  w.   We  can  therefore  write 
the  averaged  energy  spectrum  £  as  £  ■  E^^/w.   A  plot  of  £  vs 
frequency  would  give  a  smoothed  curve,  lacking  the. detail  of  the 
true  curve,  but  sufficient  to  calculate  the  total  energy  present 
over  bandwidths  comparable  to  or  larger  than  the  bandwidths  selected 
to  obtain  £"  . 


Analogously,  it  often  happens  that  what  is  presented  as  a 
pressure  spectrum  level  is  in  reality  a  smoothed  curve,  the  result 
of  plotting  data  obtained  with  filters  having  bandwidths  appreciably 
larger  than  1  Hz : 


Averaged 
spectrum 
level 

True 
spectrum 
level 


Frequency 


Use  of  too  broad  a  bandwidth  in  calculating  averaged  spectrum 
levels  will  wash  out  details.   The  spikes  in  the  true  spectrum 
level  in  the  above  example  are  lost  to  view  because  the  bandwidth 
over  which  the  averaged  levels  are  taken  is  so  wide  that  the  energy 
from  the  spike  is  small  compared  to  all  the  energy  from  the 
background  in  each  bandwidth. 
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THREE   WAYS   OF   DOING    IT 


Ipectrum  Le've^f 
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The  first  graph  presents  the  true 
spectrum  level  and  the  other  three 
present  the  same  spectrum,  except 
that  the  spectrum  levels  have 

been  averaged  over  the  bandwidth 
w.   Graph  2  presents  the 
averaged  spectrum  level  in  a 
bar  graph  format.   Graphs 
3  and  h   connect  the  data 
with  either  straight-line 
segments  or  a  smooth 
4^.    curve.   Notice  that 
Graph  k   gives  no 
information  relating  to 
the  bandwidth  of  the  filter 
obtain  the  averaged  spectrum 
this  would  have  to  be 
indicated  separately,  and  is  not 
likely  to  be  done.   Graph  3  is 
somewhat  better,  but  the  use 
of  the  straight- line  segments 
seduces  the  eye  into  assuming 
that  the  true  spectrum  level 
closely  follows  the 
straight-line  segments. 
Graph  2  is  the  least 
subject  to  being 
misinterpreted:  The 
width  of  the  individual 
bars  gives  the  bandwidth 
of  the  filter  quite 
visually,  and  each  bar  strongly 

emphasises  by  its  flat  top  that  the 
value  of  the  PSL  appropriate  to 
that  frequency  interval  is  the 
averaged  value. 

Because  of  the  fact  that 
the  averaged  spectrum  level 
smooths  out  and  erases 
detail  which  would 
otherwise  be  found  in 
-*>   the  true  spectrum  level, 
it  is  very  important 
that  the  extent  of  averaging 
over  frequency  be  well- 
This  is  accomplished  quite 
unambiguously  by  the  format 
Graph  2  or  Graph  3-   Further, 
clearly  indicates  the 

of  data  by  the  flat-topped 


Frequency 


Frequency 


indicated. 

easily  and 

of  either 

Graph  2 

averaging 


nature  of  the  curve,  whereas 
Graph  3  can  tend  to  mislead 

the  eye  in  a  casual  examination. 


4  -  „   4. 

Frequency 


SPECTRUM  LEVEL  AND  BAND  LEVEL 


Recall  that  the  total  energy  in  a  bandwidth  w  was  the  area 
under  the  energy  density  vs  frequency  curve.   The  analogous 
quantity  in  terms  of  pressure  levels  is  the  band  level  BL.   This 
is  the  expression  in  dB  terms  of  the  total  effective  pressure 
generating  the  acoustical  energy  in  the  specific  bandwidth. 

NOTICE  THAT  THE  BAND  LEVEL  IS  NOT  THE  AREA  UNDER  THE 
SPL  vs  FREQUENCY  CURVE. 

This  is  because  the  log  of  a  product  is  the  sum  of  the  logs  of 
the  individual  factors. 

Consider  a  plot  of  energy  density  £  vs  frequency  as  shown* 

£ 


It  is  clear  that  the  total  energy  E  in  the  bandwidth  w  is  given  by 

E  =  eow 
If  we  express  this  in  log  form,  we  have 


Zy 


E>,f 


+ 


It  turns  out  in  acoustics  that  taking  10- log(E/Eref )  is  the  same 
as  taking  20- log(P/Pref ),   This  means  that  the      above  equation 
is  the  same  as 

or,  treating  the  quantity  1  Hz  as  understood  but  suppressed, 


BL  =  PSL  +10  log  w 

In  general  this  expression  is  exact  only  if  the  spectrum 
level  is  constant  over  the  bandwidth  w.   However,  it  is 
reasonably  accurate  if  the  true  spectrum  level  does  not 
fluctuate  by  more  than  a  couple  of  dB  over  the  range  of  w. 


Notice  that  this  formula  can  be  turned  around t   If  a  band 
level  BL  has  been  experimentally  obtained  by  measuring  the 
acoustical  energy  over  some  bandwidth  w,  then  the  averaged 
spectrum  level  over  this  bandwidth  can  be  found  from 

PSL(averaged)  =  BL  -  10  log  w 
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COMPARISON  OF  A  SPECTRUM  LEVEL  AND  A  BAND  LEVEL  CURVE 
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Because  of  the  (10  log  w)  term  which  must  be  added  to  the 
PSL  to  obtain  the  appropriate  BL  for  a  certain  frequency  interval, 
it  is  seen  that  the  shape  of  the  BL  curve  is  such  that  the  BL 
is  stronger  at  higher  frequencies.   This  is  because  the  BL  measures 
the  total  acoustical  energy  in  some  frequency  interval,  whereas  the 
PSL  measures  the  energy  per  unit  Hz. 

Both  kinds  of  curves  have  their  advantages  and  disadvantages. 
A  curve  of  the  true  spectrum  level  gives  the  maximum  possible 
amount  of  information,  but  requires  considerable  mathematical 
analysis  if  it  is  desired  to  obtain  the  band  level  within  some 
specified  frequency  interval  from  the  spectrum  level.   A  curve 
of  averaged  spectrum  level  is  not  quite  so  difficult  (timeconsuming) 
to  use,  and  band  levels  from  it  can  be  found  more  easily,  but  often 
important  details  are  lost.   A  curve  of  band  level  is  quite  direct 
to  use,  as  long  as  the  bandwidth  agrees  with  the  particular  receiving 
system  being  used,  but  it  also  has  smoothed  out  much  detail. 
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A  SAMPLE  CALCULATION 

As  a  simple  example,  let  us  obtain  the  overall  band  level 
for  the  PSL  curve  presented  below: 

PSL 
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From  the  steplike  nature  of  the  curve,  we  can  specify  two  bandwidths 
w^  and  W2  as  shown,  with  the  associated  averaged  spectrum  levels. 

For  each  of  these,  we  can  calculate  the  appropriate  band  level* 

BL^  )  ■  B^)  +  10  log  wj 

BL(w2)  ■  B(w2)  +  10  log  W£ 

So  far,  we  have  had  to  do  no  more  than  a  little  curve-fitting  and 
application  of  formulas  which  have  been  presented  previously. 
Now,  however,  a  new  problem  has  appeared.   We  must  combine  the 
two  band  levels  into  a  single  bandlevel  for  the  total  frequency 
interval  w^  +  w2.   The  combined  BL  is  not  the  sum  of  the  individual 
ones,  because  we  are  in  dB:  Energies  add  but  logs  of  energies 
do  not.   What  must  be  done  is  to  determine  the  amount  of  energy 
given  by  BL(w^ )  and  the  amount  of  energy  given  by  BL(w2),  and 
these  energies  added  together  and  then  expressed  as  a  band  level. 

Let  us  illustrate  this  with  an  example.   We  shall  first  do 
the  example  in  the  method  indicated  above,  and  then  present  a 
simple  nomogram  which  avoids  all  the  work. 
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NOMOGRAM  FOR  COMBINING   LEVELS 
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This  nomogram  allows  the  combination  of 
two  different  band  levels  (or  two 
independent  sound  pressure  levels)  by 
graphical  means.   If  the  two  levels  to 
be  combined  are  L^  and  L£,  the  procedure 
is  as  follows: 

Obtain  the  difference  L^  -  L2- 

Locate  this  number  on  the  horizontal 
axis,  move  straight  up  until  the 
— I  curve  is  reached,  and  then  move 
horizontally  to  the  left  until  the 
vertical  axis  is  reached,  and  read 
that  number. 


Add  that  number  to  L^ 
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As  an  example  of  the  application  of  th 
problem  worked  out  on  the  previous  pag 
which  were  calculated  there  were  BL^  ■ 
BLo  =  168  dB  re  l^Pa.   The  difference 
value,  projecting  upward,  and  reading 
so  that  the  combination  of  these  two  1 
band  level  of  170  +  2.1  =  172.1  dB  re 
uncertainty  in  levels  exceeds  a  few  te 
appropriate  to  round  this  combined  lev 


is  nomogram,  recall  the 
e.   The  two  band  levels 
170  dB  re  l/"Pa  and 
is  +2  dBj  finding  this 
to  the  left  gives  2.1  dB, 
evels  yields  a  total 
l^Pa.   Since  the  usual 
nths  of  a  dB,  it  is 
el  to  171   dB  re  l^?a. 


TONALS    COMBINED  WITH  BACKGROUND 
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When  there  are  strong  tonals  present,  the  band  level  cannot 
"be  expressed  as  B  +  10  log  w,  because  of  the  additional  energy 
contributed  by  the  spike(s).   We  must  calculate  the  energy  from 
the  background  spectrum  (without  the  spikes)  and  add  this  to  the 
energy  of  the  monofrequency  spikes  which  are  also  present: 


BLi 


B  +  10  log(w-l ) 


+  10  log  w   if  w  2>  1 


BLA  =  A  +  10  log  1 


The  total  band  level  is  found  by  combining  BL^  and  BLg  with  the 
help  of  the  nomogram. 

Notice  that  if  the  tonal  is  relatively  strong  and  the  bandwidth 
narrow  enough,  the  combined  level  will  be  due  essentially  to  the 
tonal.   As  the  bandwidth  w  increases,  however,  the  influence  of 
the  (10  log  w)  term  will  become  stronger,  causing  BLg  to  increase. 
Eventually  the  contribution  from  the  background  will  get  sufficiently 
larsre  that  BL^  will  be  sufficiently  large  that  the  influence  of 
BL^  will  be  lost. 
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T  HR  SONAR  EQUATION 


SL  >ML  +  DT 


EL  =  Echo  Level  for  active  sonar  or  the  received  Signal  Level  for  passive 
sonar.  It  is  the  level  in  dB  of  the  desired  sifnal  vhen  it  reaches 
the  receiver. 

ML  ■  Masking  Level.  It  is  the  level  in  dB  of  the  undesired  signal  at  the 
receiver  -which  is  competing  vith  the  Echo  Level. 

DT  ■  Detection  Threshold.  This  is  the  number  of  dE  by  vhich  the  Echo  Level 
must  exceed  the  Masking  Level  to  allov  the  receiving  system  to 
register  a  detection  with  a  specified  degree  of  confidence.  This 
quantity  involves  the  probability  of  detecting  the  target  in  the 
presence  of  noise  and  also  the  probability  of  giving  a  false  alarm. 


The  SOIj'AR  equation  states  that  signal  at  the  receiver  must  equal  or  exceed 
backf round  noise  by  an  amount  specified  by  the  Detection  Threshold  if  a 
specified  decree  of  performance  is  to  be  achieved. 
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LIST  CF  SYMBOLS 


EL  • 

•  Echo  Level 

ML  « 

■  Masking  Level 

DT  ■ 

•  Detection  Threshold 

SL  « 

■  Source  Level 

TL  • 

■  Transmission  Loss 

PSL  « 

«  Pressure  Spectrum  Level 

BL  - 

»  Band  Level 

SSL  • 

«  Source  Spectrum  Level 

AG  • 

»  Array  Gain 

NL  ■ 

>  Noise  Level 

NSL  « 

»  Noise  Spectrum  Level 

DI  « 

■  Directivity  Index 

w   «■  bandwidth 

d   -  detection  index 

wt  ■  time -bandwidth  product 
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PASSIVE  SONAR 


r~T) 


L 


6 


EL     =     SL  -  TL 


SL     -     TL     X    ML     +     DT 


>' 


where     SL  ■     Source  Level  of  the  target  . 
TL  =  One-way  transmission  loss 

xiach  term  in  the  above  equation  will  now  be  discussed  separately. 
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SLJ-  TL  ^  ML  +  DT 


SOURCE  LEVEL  -  PASSIVE  SONAR 


Tonals  (monfrecuency  components) 

Sources:   Propulsion  machinery  (Main  motors,  reduction  [ear) 

Auxiliary  machinery  (generators,  pumps,  air-conditioners) 
Propeller-induced  resonant  hull  excitation 
j iydrodynamically -induced  resonant  excitation  of  cavities, 
plates,  and  appendages 


SL  =  Level  of  acoustic  ener£y  1  meter  from  tie  acoustic  center  of  the 
source. 

measurements  are  usually  made  at  a  distance  from  the  source  and 
extraoolated  back  to  1  meter. 


II .  Proadband  noise 

Sources:   Cavitation  at  or  near  the  propeller 
Radiated  flow  noise 
Caviation  at  struts  and  appenda£es 


SL  -  SSL  +  10  lo£  w 

■where  "w  is  the  bandwidth  of  the  system  and 

SSL  =  pressure  spectrum  level  1  meter  from  the  acoustic  center 
of  the  source. 

If  the  SSL  is  not  essentially  constant  over  "w,  the  bandwidth 
must  be  subdivided  into  smaller  segments,  the  PL  of  each 
found,  and  then  all  the  PL's  combined  to  ^ive  SL 
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TYPICAL  SOURCE  SPECTRUM  LEVELS 


Typical  source  spectrum  levels  for  a  conventional  submarine  on  electric 
drive  at  periscope  depth.  SSL  re  If Pa. 


freq  (Hz) 
Speed  (knots, 


k 
6 
8 
10 


100 


1000 


133 
LUO 

LU3 
11*7 


120 
127 
133 
135 


10,000 


100 
112 

119 
122 


Adapted  from  Urick,  Principles  of  Underwater  Sound  for  Engineers 
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SL  -  TL    >(ML)+  DT 
MASKING   LEVEL 


?  \ 


4 


Unwanted  noise  may  come  in  uniformily  from  all  direction  or  it  may  arrive 
from  a  pre fere d  direction. 


where 


ML  -  NL  -  AG 


ML  »  masking  iciel 

NL  ■  noise  level  in  the  environment 

AG  ■  array  gain,  the  improvement  (reduction)  in  dB  against  the 
noise  field  resulting  from  the  array  design. 


70 


ML  -  NLj  -  AG 
Noise  Level 


! 

NL  =  NSL  +  10  lot  w 


where 


NL  •  noise  level,  the  acoustic  energy  in  dP  found  in  the  bandwidth 

of  the  system  , 

NSL  ■  noise  spectrum  level,  pressure  spectrum  level  of  the  unwanted 

noi  se , 

and 

w  *  the  bandwidth  of  the  system  in  Hz. 


Kote:   The  same  limitations  apply  in  using  this  formula  as  were 
pointed  out  in  the  discussion  of  the  broadband  SL. 
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NL 


(NSl)  +  10  log  w 


NOISE  SPECTRUM  LEVEL  (Deep  water) 


100 


100 
Frequency  (Hz) 


Adapted  from:  Wenz,  Jour.  Acoust.  Soc.  Am.  3k   ,  1936-1956  (1962) 
Perrone,  ibid  1*6,  762-770  (19697 
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ML  =  (MSL  +  10  lo[  vj 


N0IS3  SP3CTKDM  L3V3L  (Self-noise) 


Sources  of  self -noise:   J iydro dynamic  noise  due  to  flov  past  receiver 

Acceleration  of  receiver 
Machinery  noise  of  ov:n  "ship" 


Values  usually  presented  as  an  "equivalent  omnidirectional"  noise  spectrum 
level. 


KL( self -noise)  »  KSL( self -noise)  +  10  loj.  w 


If  significant  self -noise  is  present,  it  must  te  combined  with  the  ambient 
noise  level  to  obtain  the  total  noise  level.  Since  these  two  sources  of 
noise  are  independent,  their  energies  add  and  the  nomogram  can  be  used. 
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ML   =   NL  -.'AG 


ARRAY  GAIN  AND  DIRECTIVITY  INDEX 


Array  gain  is  the  decrease  in  masking  level  from  the  noise 
level  brought  about  by  the  design  of  the  receiving  array. 


If  the  noise  is  isotropic  (arriving  uniformly  from  all  directions) 
and  the  signal  arrives  from  one  direction  only,  then 

AG  =  DI 

where 

DI  =  Directivity  Index 


=  10  log  (the  ratio  of  the  noise  energy  sensed 

by  an  omnidirectional  receiver  to  that 
sensed  by  the  directional  receiver) 


In  this  case, 


ML  =  NL  -  DI 


Array  gain  can  include  both  the  design  of  the  receiving  array 
and  also  the  signal  processing.   Signal  processing  such  as  correlation, 
by  taking  into  account  the  different  spatial  or  temporal  properties 
of  the  signal  and  noise  can  result  in  an  array  gain  different  from 
the  directivity  index. 
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4 

4- 

Idd  together,  + 
jet  N  times      4.  ^ 
;he  signal 

+ 


■  suiting  response 
is  a  function  of 
ingle: 


Pair  by  pair, 
the  signals 
cancel 


*& 


'olar  presentations 


weak 
reception 


.strong 
reception 


20  log  (P/P— ^ 


-  a* 


zyo 


»I  denotes  how  directional  the  receiver  is. 

DI  -  0  dB     Omnidirectional  receiver 

DI  ■  20  dB     ReceiTing  only  1%  of  the  energy  present  in  the  noise 
field. 

1  depends  on  sise  of  the  array  compared  to  the  wavelength.  For  example, 
or  a  long  line  array  of  length  L  much  greater  than  ^ 

DI  -  10  log  (2L/<\  ) 

•nerally,  as  frequency  increases  the  receiver  becomes  more  directive  (DI  Increases) 
lit,  if  the  frequency  is  too  high  we  may  get  more  major  lobes. 
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A  PHYSICAL  MEANING  OF  DIRECTIVITY  INDEX 


Sphere  of 
unit  radius 


Illuminated 
patch  on  the 
unit  sphere 


Surface  area  of  the  unit  sphere  *  ^TT(l) 

Area  of  the  illuminated  patch  ■  _fV.(0 

where  St   is  the  solid  angle  subtended  by  the  patch. 

The  fractional  area  of  the  sphere  which  is  illuminated 
is  given  by  SL  /kV   and,  for  this  simple  example,  the  directivity 
D  is  given  by  u 

The  directivity  index  is  the  expression  of  the  directivity  D  in 
dB  terms,       r 


DI  =  10  JZo*  D 


A  more  general  definition  of  the  directivity  is  given  as 
follows:   If  the  pressure  field  generated  by  an  array  acting  like 
a  source  (or  the  relative  sensitivity  of  an  array  acting  like  a 
receiver)  is  written  in  the  form 


P(r,O,0) 


Pax(r)  H(O,0) 


where  Pax  is  the  pressure  on  the  acoustical  axis  of  the  array, 
then  H(9,0O  is  normalized  to  have  a  maximum  value  of  unity  on 
the  acoustical  axis,  and  to  be  no  greater  than  unity  for  any  other 
direction.   The  directivity  can  then  be  expressed  as 
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SOME  ROUGH  ESTIMATIONS  OF  DI 


If  we  have  a  line  array  with  a  single  major  lobe,  and  if 
the  angular  width  of  the  major  lobe  is  0  (in  radians),  then  if 
the  array  is  surrounded  by  a  unit  sphere,  the  illuminated  patch 
appears  as  below 


and  the  directivity  can  be  crudely  estimated  to  be 


(  0  in  radians  ) 
where  it  must  be  understood  that  0  is  small. 


A  little  more  careful  estimation  of  the  directivity  from 
the  geometry  of  the  radiation  (or  reception)  pattern  leads  to 
some  relatively  simple  but  only  approximate  estimates  for  the 
directivities  of  arrays  of  simple  design, 

Line- like  arrays: 


D  ~  f  tzxj&ih    »  , 


Piston-like  arrays: 


Z~{Tf j  J  »! 


These   formulas  can  be  used  as   rough  estimates   if  the  arrays  are 
highly  directive   (DI   in  excess  of  about  6  dB  at  the   least),    are 
designed  to  have  only  one  major  lobe,    and  are  not   strongly  shaded 
or  steered. 

An  array  is  shaded  when  different  elements  have  different  amplifications. 
This  is  done  to  alter  the  width  of  the  major  lobe  and  the  relative  strengths 
of  the  side  lobes.     Narrowing  the  major  lobe  will  increase  the  strengths  of 
the  side  lobes,   and  vice  versa.     An  array  is  steered  to  point  the  major  lobe 
in  various  directions  with  respect  to  the  array.     This  is  done  by  introducing 
time  delays  in  the  outputs  of  the  elements. 
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PASSIVE  SONAR  EQUATIONS 


Let's  put  it  all  together: 


BROADBAND  WITH  NEGLIGIBLE  TONALS 


SSL  ♦  10  log  v     -     TL   >     NSL    *     10  log  w     -  AG     «■     DT 


— 

SSL 

-  TL 

> 

NSL  • 

-  AG 

♦ 

DT 

TONALS  WITH  NEGLIGIBLE  CONTRIBUTION  FROM  BACKGROUND 


SL  -  TL   >  NSL    +     10  log  w     -     AG     ♦     DT 
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FIGURE  OF  MERIT 
(FOM) 

The  passive  SONAR  equations  can  be  rewritten  by  solving  forTL: 


where 


or 


TL      <C     FOM 


FOM  -  SL  -  (NSL  +  10  log  w  -  AG  +  DT) 


FOM  -  SSL  -  (NSL  -  AG  +  DT) 
for  totoal  or  broadband  detection  respectively. 
Graphically: 


TL  V 


FOM  " 


A 


range 


ranges  for  which  detections  are  possible  at  conditions 
equal  to  or  better  than  specified  by  DT. 
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To  take  advantage  of  a  narrow  bandwidth  and  still  permit  reception  of  a 
Doppler  shifted  signal,  it  is  necessary  to  design  a  system  that  employs 
many  parallel  filters  which, while  individually  narrow, cover  the  entire 
frequency  range  of  interest. 

The  output  of  such  a  system  might  be  presented  in  this  manner: 


a 

1 

<D 

c — L^J ' 1 

i   i 1 

tl  t2 


Time 


At  any  instant,  the  output  of  each  of  the  individual  filters  might  look 
like  this: 


At  t,: 


At  t2: 


-p 

s. 
1 


T---^ 



I_r 

frequency 


I 

o 


LP-dl 


frequency 


Att.     | 

o 


frequency 


where  T  is  the  threshold  setting. 
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SL  -  TL  <  ML  +  @ 
DETECTION  THRESHOLD 


The  detection  process  consists  of  designating  a  threshold  which,  when  exceed, 
causes  a  detection  to  be  recorded.  If  the  signal  is  much  stronger  than  the 
noise,  it  is  clear  that  a  threshold  can  be  defined  that  will  allow     valid 
signals  to  be  recorded  while         r  ignoring  the  noise.  However,  when  the 
signal  and  noise  are  of  comparable  size,  any  threshold  that  will  catch  a 
reasonable  number  of  valid  signals  will  also  record  "detections"  when  a  valid 
signal  is  absent. 


The 
for 


probablity  of  obtaining  a  certain  amplitude  of  the  noise  alone  as  well  as 
the  signal- and- noise  is  sketched  below: 


\  detection 


'//y    false  alarm 


s&n 


amplitude  A 


where 

Ajj  is  the  mean  value  of  the  amplitude  with  noise  alone 
A  s&n  is  the  mean   value  of  the  amplitude  with  signal  and  noise 
<r  is  a  measure  of  the  width  of  the  curves  (assumed  equal) 
T  is  the  threshold 

The  shaded  area  marked  "detection"  is  the  probablity  of  making  a  valid  detection 
when  the  threshold  is  T.  (The  area  under  the  curve  is  taken  as  unity.)  And 
the  shaded  area  marked  "false  alarm"  is  the  probablity  of  making  a  "detection 
when  the  signal  is  absent* 

Let  P(D)  ■  the  probability  of  detection 

P(FA)  -  the  probability  of  a  false  alarm 
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DETECTION  INDEX 


(D 


Amplitude 


^s+n  ~  ^n 


h 


S1« 


We  need  to  find  some  simple 
criterion  which  will  describe 
how  the  curves  of  noise  alone  and 
signal  and  noise  combined  overlap, 
since  for  a  given  value  of 
threshold,  it  is  the  overlap  of 
the  curves  which  establishes 
the  relative  values  of  P(D)  and 
P(FA). 

Notice  the  configuration  of  the 
two  curves  given  in  the  first 
graph.  The  curves  are  quite  well- 
separated,  in  that  their  peaks 
have  amplitudes  whose  difference 
is  larger  than  the  o*  .  These 
two  curves  could  be  considered 
well-eeparated . 

In  the  second  graph,  each  curve 
has  the  same  average  value  of 
A  as  before,  but  now  the  0*' is 
twice  as  large.  The  curves  now 
are  not  well- separated. 
In  the  third  graph,  the  curves 
have  the  same  0*'  as  in  the  second, 
but  the  two  have  been  shifted 
further  apart  so  that  they  are 
well-separated.  Indeed,  compariso 
of  the  first  and  third  graphs 
will  reveal  that  the  amount  of 
overlap  of  the  pairs  of  curves 
is  the  same.  Notice  that  the 
ratio  of  the  difference  in  the 
iverage  values  of  the  amplitudes 

first  graph  is  identical  with 
that  for  the  third.  Thus,  it 
is  the  separation  between  the 
two  curves  ih  units  of  the  width 
of  the  curves  which  measures 
how  well- separated  the  two  curves 
are.  It  is  thus  plausible  to 
define  the  detection  index  d  by 


</  = 


A***  "*  A  * 
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s  +  ir\ 


s+io 


H* 


Let  us  study  the  effects  of  changing  the  positions  of  the  threshold  and  the 
detection  index  on  the  probabilities  of  detection  and  faDLse  alarm.  If  the 
original  configuration  of  the  noise  curve,  signal  and  noise  curve,  and 
threshold  are  as  given  in  the  top  curve,  then  we  can  establish  the  indicated 
P(D)  and  P(FA)  by  the  shaded  areas.  Now,  if  the  threshold  is  lowered,  as  in 
the  second  figure,  then  P(D)  shows  some  increase,  but  the  P(FA)  increases 
rather  drastically:  lowering  the  threshold  will  increase  the  probability  of 
a  detection,  but  the  "clutter"  will  also  increase.  On  the  other  hand,  if  the 
noise  curve  and  threshold  are  maintained  constant,  but  the  signal  gets  stronger, 
then  the  curve  of  signal  and  noise  will  move  to  the  right.  This  means  that  the 
detection  index  increases,  and  the  P(D)  increases  while  P(FA)  remains  the  same. 
This  says  no  more  than  that  a  stronger  signal  has  a  greater  probability  of 
being  detected. 

If  the  statistical  properties  of  the  noise  and  the  combined  signal  and  noise 
are  known,  then  it  is  possible  to  make  curves  of  the  P(D)  and  P(FA)  relationships 
for  various  values  of  detection  index.  These  are  known  as  Receiver--'  Operating 
Characteristics  (ROC)  curves. 
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RECEIVER  OPERATING  CHARACTERISITCS 
ROC 


1*0 


Q 

"5T 


0.2        O.li 


0.6 


0.8 


1.0 


P(FA) 


Typical  ROC  Curves 
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FALSE  ALARM  RATE 

The  rate  of  false  alarms  is  the  product  of  the  sampling  rate  and  P(FA). 
Since  sampling  rates  are  generally  high,  the  importance  of  a  very  low 
P(FA)  should  be  noted. 

For  example:  If  the  output  of  each  of  10  parallel  filters  is  sampled 
for  one  second  each,  a  P(FA)  of  l£  will  result  in  an  average  of  one 
false  alarm  every  10  seconds. 

The  sketch  below  shows  a  representative  plot  of  the  relationship  between 
d  and  P(FA),  for  a  P(D)  -  0.5; 


30  r 


P(D)  -  0.5 


P(FA) 
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ENERGY  DETECTION 

For  a  passive  detection  system, 

a  completely  unknown  signal,  and  Gaussian  noise,  best  performance  is 
obtained  by  filtering  in  a  bandwidth  w,  taking  the  square  of  the  signal, 
and  accumulating  over  a  processing  time  t.  If  the  signal-to-noise 
ratio  is  small,  it  can  be  shown  that  the  detection  threshold  for  this 
process  is 

DTe  -  $   lot  (d/wt) 

where  wt  is  called  the  time-bandwidth  product. 

To  decrease  DT  for  a  fixed  d,  it  is  necessary  to  increase  the  time -bandwidth 
product. 


Footnote : 

The  detection  threshold  used  in  this  material  is  that  calculated  over 
the  complete  bandwidth  of  the  processor.  It  often  occurs  that  what  is 
specified  is  the  detection  threshold  in  a  one  Hertz  bandwidth,  DT(1). 
This  latter  quantity  occasionally  leads  to  SONAR  equations  which  have  a 
slightly  simpler  format  when  it  is  assumed  that  the  spectrum  levels  are 
constant  over  the  bandwidth  of  the  system.  Otherwise,  we  feel  that  the 
use  of  DT(1)  introduces  more  possibilities  for  confusion,  so  that  its  use 
creates  more  problems  than  it  solves.  (In  particular,  if  Urick's  excellent 
book  is  studied,  it  must  be  fully  understood  that  whereas  he  uses  DT  throughout 
most  of  the  text,  his  chapter  devoted  to  the  detection  threshold  uses  DT(i).) 
The  formula  relating  the  two  quantities  is 
• 

DT(1)  -  DT  +  10  log  w 
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DOPPLER   SHIFT 
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While  decreasing  w  will  improve  the  ability  of  a  paasive  SONAR  to  detect 
tonals,  it  cannot  be  made  smaller  than  the  expected  Doppler  shifts. 
For  a  passive  SONAR  system ,  this  requires  that  the  minimum  acceptable 
bandwidth  be 

w(in  Hz)  =  0.7  x  (frequency  in  kHz)  x  (range  rate  in  knots) 
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EXAMPLE  1 

A  submarine  with  a  160  dB  re  1 ^  Pa  line  at  $00  Hz  crosses  a  convergence 
zone  at  kO   km  from  an  omnidirectional  hydrophone.  The  bandwidth  of  the 
receiver  is  100  Hz  and  the  sea  state  is  3»  Assuming  a  convergence  gain 
of  12  dB,  find  out  how  long  the  signal  must  be  accumulated  to  give  a 
P(D)  -  0.50  and  a  P(FA)  -  10"**. 

Solution: 

TL  =  20  lot  r     ♦     ar     -  G     =     20  log  U  x  10^)  ♦  $  x  10"*   •  h  X  1&     -     12     -     82  dE 

ML  -  NSL  +  10  log  w  -  DI     -     66+20-0     «     86  dB 
SL  -  TL   >   ML  «•  DT 

DT  .v.   SL  -  TL  -  ML    =  160  -  82  -  86    -     -8  dB 
DT  -  5  log  (d/wt) 
where       d  •  15 

log  (wt/d)  -  1.6 
wt/d  -  kO 
wt  -  600 
t  ■  6  sec 
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EXAMPLE  2 


Sams  as  Example  1,  but  let     w    be  just  vide  enough  to  encompass  a  20  knt 
target. 

Solution: 

w  ■  0.7  x  (frequency  in  kHz)  x  (range  rate  in  knots) 

-  0.7  h    20  -  •,  7  Hz 

ML  »  NSL  ♦  10  log  v  -  DI  -   66  +  8.5  -  71*. 5  dB 

DT  <  SL  -  TL  -  ML  -  160  -  82  -  71*. 5  -  3.5  dB 

DT  -  5  log  (d/vt) 

where   d  -  15 

log  (d/wt)  =  0.7 

d/wt  -  5 

wt  •  3 

t  -  3/7  sec 
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EXAMPLE  3 

A  conventional  submarine  at  periscope  depth  is  traveling  at  4  kts.   The 
Sea  State  is  3  and  the  layer  depth  is  100  m.   At  what  range  will  the 
submarine  be  detected  by  a  36  m  deep  hydrophone  listening  at  1000  Hz 
if  DI  =  20  dB  and  DT  =  0  dB? 

Solution: 

SSL   -   TL  >,      NSL  -  DI  +   DT 

TL  <T   SSL  -   NSL  +  DI  -   DT  =   120  -62  +  20   =   78  dB 

For  sound  to  be  trapped  in  the  duct, 

f  >   f  . 

'         min 


5   3/2 
=   2-10  A>     =   200  Hz 


So  we  use  the  TL  for  a  mixed  layer 

TL  =  10  log  r^  +  10  log  r  +  ar  +  b (r/r  ) 
t  s 

r   =   1.05-102  D/(d  -  z  )  *' 2  =   1 .05 "102 (100/8)   =   1310  m 
t  s 

2   1/2  2 

r   =   8.4-10  D  '         =   8.4-10  -10  =   8400  m 

s 

From  the  graphs  and  equations 

-5 
a  =  6-10    dB/m 

b  =  1.7  dB /bounce 

-4 

Notice  that   ar  +  b(r/r  )   +   (a  +  b/r  )  r  =   2.6*10   r 

s  s 


Thus     78  =  10  log  1310  +  10  log  r  +  2.6-10 

10  log  r  +   2.6-10~   r  =   47 
Now,  trial  and  error 


-4 


r 

10  log  r 

2.6-10"4  r 

TL 

5-104 

47 

13 

60 

2-104 

43 

5.2 

48 

4 
1.5-10 

41.8 

3.9 

46 

4 
1.7-10 

42.3 

4.4 

47 

r  =   1.7-10  m 
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EXAMPLE  h 


Repeat  Example  3  but  with  DI  «  0  dB. 


Solution: 


TL   £   SSL     -     NSL     ♦     DI     -     DT     -     120  -62     •     $8  dB 

-4 
TL     -     10  log  r       ♦     16  log  r     +  2.6  x  10       r 

58  «  10  log  1310  +  10  log  r  +   2.6  x  10-4  r 

58  »  31  ♦  10  log  r  +   2.6  x  10"4  r 

10  log  r  +  2.6  x  10~4  r  =   27 


10  log  r   2.6  x  10 


-4 


5  x  io2 


27 


0.13 


TL 


27 


Bingo 


But    r  <  rt 

So  we  must  use  this  form  of  the  TL; 

TL  »  20  log  r   ♦  ar 

58  -  20  log  r 

log  r  -  2.9 

r  -  800  m 
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EXAMPLE  $ 


Repeat  Example  3  with     f  •  100  Hz.       Assume  "heavy  shipping" as  would  be 
encountered  in  the  Mediterranean  Sea. 


Solution: 


f    <  f 


min 


No  trapping  in  the  layer.     Use 

TL     -     20  log  r     +  ar 

TL    ^     SSL     -     NSL     +     DI 

20  log  r    +     10"6  r    -     77 


-     DT     -     133-76    +20    -     77  dB 


r 

20  log  r 

10"6  r 

TL 

8  x  103 

78 

8  x  10*3 

78 

7  x  103 

77 

7  x  10'3 

77 

7  x 


103  m 


Bingo 
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ACTIVE  SONAR 


SL  -  TL 


SL  -»  TL  +  TS 


-  TL  ♦  TS  -  TL 


SL 


2TL  +  TS  ^  ML  ♦  DT 


where 


SL  -  source  level  of  transmitter 
2TL  •  two-way  transmission  loss 
TS  m     target  strength 


Each  term  in  the  above  equation  will  now  be  discussed  separately. 

Note:  the  above  active  SONAR  equation  is  for  a  monostatic  case, 
wherein  the  source  and  receiver  are  coincident  so  that  the  TL 
from  source  to  target  is  the  same  as  the  TL  from  target  to 
receiver.   Our  discussion  will  usually  be  restricted  to  this 
case.   Extension  to  the  bistatic  situation  in  which  source  and 
receiver  are  at  different  locations  is  trivial:   If  the 
transmission  loss  from  source  to  target  is  TL,  and  from  target 
to  receiver  TL' ,  then  all  that  is  required  is  to  replace 
2TL  with  TL  +  TL' ,  so  that  we  have 


SL 


TL   -   TL' 


TS 


> 


ML 


DT 
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LIST  OF  ADDITIONAL  SYMBOLS  FOR  ACTIVE  SONAR 

SL  »=  source  level  of  the  active  sonar  system 

TS  ■  target  strength  of  the  desired  target 

RL  ■  reverberation  level 

TS'  ■  target  strength  of  the  scatterers  causing  reverberation 

TL  *  transmission  loss  from  source  to  target 

TL'  =  transmission  loss  from  target  to  receiver,  equal  to  TL  in 
the  case  of  monostatic  systems 


Le  =  "effective  length"  of  a  target 

t  ■  pulse  duration  of  the  acoustical  signal  from  the  source 

sv  -  a  measure  of  the  scattering  cross-section  per  unit  volume 

sg  *  a  measure  of  the  scattering  cross-section  per  unit  area 

Sv  *  scattering  strength  for  volume  scatterers 

S_  ■  scattering  strength  for  surface  scatterers- 
s 
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SL)  -  2TL  +  TS  \   ML  +  DT 


7 


SOURCE  LEVEL  -  ACTIVE  SONAR 


SL  »  sound  pressure  level  1  m  from  the  acoustic  center  of  the 
transmitting  transducer. 


SL  -  16  log  (Pover)  +  DI  +  171 


where 


or 

where 
and 


Power  ■  Acoustic  power  output  of  the  transducer  in  Watts 

DI  =  10  log  (the  ratio  of  the  intensity  on  the 
acoustic  axis  to  the  intensity  at  the  same 
distance  from  an  omnidirectional  transmitter 
with  the  same  acoustic  power  output) . 


DI  -  20  log  P(r)/P0(r) 

?oM       m      pressure  produced  by  an  omnidirectional  transducer 
with  the  same  power  output 


171     ■    factor  to  take  care  of  the  units 


u 
« 

a> 

© 
o 

I 

CO 


2U0  r 


220 


200 


180 


160 


- 

i  a                 ^" 

- 

K«   30  « 

^^ 

^^~~^ 

\*^*~^ 
<>^~^^ 

c^— 

0.1 


1         10         100      1000 
Acoustic  power  output  (Watts) 


10,000 
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SL  -  2TL  +  (TS)^  ML  +  DT 
TARGET  STRENGTH 


TS  ■  (sound  pressure  level  of  reflected  wave  1  m  from  the 
acoustic  cenMr  of  the  target  in  a  direction  back 
towards  the  receiver)  -  (sound  pressure  level  of 
the  incident  wave  at  the  position  of  the  target  if 
the  target  were  absent). 

or 

TS  -  20  log  P'(r»-1)/P$r) 


where 


P'(r'»l)  is  the  pressure  of  the  reelected  wave  1  m  from  the 
acoustic  center  of  the  target  in  a  direction  back  towards  the 
receiver. 


Target  Strength  depends  on  shape,  construction,  and  orientation  of  the 
target,  and  on  the  frequency  and  pulse  length  of  the  sound  wave. 
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EFFECT   OF    TARGET    SIZE 


\^/eff<?c4iv'<?   ie»*ii,     L 


+-. -JV-^^-^/V^-AA- 


<-—  c±  — H 


Because  of  the  effective  length  L  of  the  target,  portions  of  the 
incident  pulse  which  reflect  from  the  tail  (in  the  above  sketches) 
travel  a  total  distance  of  2L   further  than  the  reflections  from 
the  nose.   This  means  that  the  apparent  duration  of  the  pulse 
is  stretched  from  t  to   (t  +  2L  /c)  .   For  very  long  pulses, 
the  additional  time  interval  corresponding  to  the  round  trip, 
2L  /c,  can  be  neglected.   For  short  pulses,  however,  this  stretching 
can  be  quite  important  in  determining  the  target  strength. 
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The  Table  shows  a  few  representative,  but  very  approximate,  values  of 
target  strength  for  a  variety  of  objects, 


Target 

Frequency 
in  kHz 

Fleet  Submarine 
(approx  100  m  long) 

10 
2k 

S  Class  Submarine 
(approx  70  m  long) 

2k 

Torpedo 

2k 

Surface  vessel 

3  ft  mine 

30-90 

Approximate  Target  Strength  (dB) 
Beam  Aspect        Bow-Stern  Aspect 
35  30 

25  16 

20  13 


10  to  30 
-8 


-20 


Extracted  from:     Physics  of  Sound  in  the  Sea,  NDRG 

Principles  of  Underwater  Sound  for  Engineers,   Urick 
Principle  and  Applications  of  Underwater  Sound,  NDRC 
NAVSHIPS  0967-129-3010 
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SL   -   2TL 


+   TS  ^  (ML)  +   DT 


MASKING  LEVEL 

The  undesired  noise  competing  with  the  desired  signal  can 
arise  from  two  possible  mechanisms  in  the  case  of  active  SONAR: 

Masking  from  ambient  or  self  noise 

Masking  from  reverberation 

Which  of  these  two  mechanisms  is  the  dominant  depends  on  the 
range  between  source  and  target.   Generally  speaking,  for  low- 
powered  active  systems  the  performance  will  be  noise-limited, 
in  that  the  maximum  effective  detection  range  will  be  determined 
when  the  echo  level  diminishes  below  that  level  for  which  it  can 
be  extracted  from  the  ambient  or  self  noise.   For  higher-powered 
active  systems,  however,  reverberation  becomes  an  important  source 
of  interference.   Reverberation  tends  to  decrease  with  increasing 
range,  but  less  rapidly  than  the  echo  level.   Thus,  if  the  echo 
level  has  diminished  until  it  is  buried  in  the  reverberation, 
but  the  reverberation  has  not  yet  fallen  sufficiently  below 
the  ambient  (or  self)  noise,  then  the  maximum  detection  range 
is  determined  by  the  reverberation.   The  two  situations  are 
suggested  in  the  figures: 


Low-Powered 
System 


range 


The  choice  is  a  little  more  subtle  than  suggested  by  the  simple 
curves  above.   The  detection  threshold  for  noise  masking  may 
be  different  than  that  for  reverberation  masking,  And  the 
detection  threshold  for  reverberation-limited  performance  may 
depend  strongly  on  the  doppler- induced  difference  between  the 
frequency  of  the  reverberant  sound  and  that  reflected  from  the 
desired  target. 
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AMBIENT  AND  SELF  NOISE 

The  masking  level  is  given  by 

ML  ■  NL  -  AG 

where  the  noise  level  NL  is  that  from  the  ambient  noise  of  the 
ocean,  noise  generated  by  the  receiving  platform,  or  the  combination 
of  the  two  if  they  are  of  nearly  equal  strengths. 

Recall  that  if  the  noise  spectrum  level  is  sufficiently 
smooth  over  the  receiver  bandwidth,  then  we  can  write 

NL  =  NSL  +10  log  w   . 

Otherwise,  the  total  bandwidth  must  be  broken  down  into  smaller 
bands,  the  band  level  found  for  each,  and  then  these  band  levels 
combined  according  to  the  nomogram  for  combining  independent 
levels. 
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REVERBERATION 

If  the  sound  field  generated  by  the  source  is  strong  enough,  and 
there  are  enough  scatterers  present  in  the  vicinity  of  the  desired 
target,  then  unwanted  reflections  and  scattering  of  sound  from 
the  undesired  scatterers  are  received  along  with  the  echo  from 
the  desired  target.   The  echo  level  of  the  signal  received  from 
the  desired  target  is 


EL  ■  SL  -  2TL  +  TS 


(monostatic) 


where  TS  is  the  target  strength  of  the  desired  target.   The 
reverberation  level  RL  can  be  written  in  the  form 


RL 


SL  -  2TL  +  TS' 


(monostatic) 


where  TS'  is  the  combined  target  strength  of  all  the  undesired 
targets  which  are  generating  the  masking  level.   Since  the  RL 
is  the  masking  level  when  reverberation  is  the  dominant  source 
of  unwanted  sound,  the  SONAR  equation  becomes 


or 


EL  -  RL  >   DT 

TS  -  TS'  ^  DT 


(monostatic) 


The  central  problem  in  obtaining  TS'  is  determining  the  volume 
or  surface  area  containing  the  scatterers  which  is  illuminated 
by  the  source,  sending  sound  to  the  receiver,  and  sending  the 
sound  so  that  it  arrives  at  the  same  time  as  the  signal  from  the 
desired  target.   This  is  suggested  in  the  sketch  below. 


Reverberation- 
Limited  Performance 
(bistatic) 
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VOLUME  REVERBERATION 


The  target  strength  TS'  of  the  undesired  scatterers  can  be 
written  as 


TS 


10  log  SVV 


where  V  is  the  volume  within  which  the  reverberation  sent  to 
the  receiver  is  in  competition  with  the  reflected  signal  from 
the  desired  target  and  sv  is  a  measure  of  the  scattering  cross 
section  per  unit  volume.   It  has  become  conventional  to  express 
sv  in  dB  terms  through  the  definition 

Sv  ■  10  log  sv 

where  Sv  is  the  scattering  strength.   Thus,  we  have 


TS' 


Sv  +  10  log  V. 


The  calculation  of  the  reverberating  volume  V  contains  the  heart 
of  the  problem.   As  suggested  previously,  this  is  a  study  in 
geometry,  involving  the  way  in  which  major  lobes  of  source  and 
receiver  overlap. 

If  we  restrict  ourselves  to  the  monostatic  situation,  and  assume 
that  the  major  lobe  of  the  receiver  completely  overlaps  the 
major  lobe  of  the  source,  it  then  follows  that  the  volume  V  is 
given  by 

10  log  V  -  10  log  pL     +  TLg  +10  log  SI 

where  t  is  the  time  duration  of  the  acoustical  pulse  and  IL    is 
the  solid  angle  subtended  by  the  major  lobe  of  the  radiation 
pattern  of  the  source.   Recall  that  this  is  related  to  the 
directivity  index  by  DI  ~  10  log(^7T/fV).   The  term  TLg  arises 
from  the  fact  that  the  cross  sectional  area  of  the  reverberating 
volume  clearly  depends  on  the  range  from  the  source,  and  this  can 
be  seen  to  be  expressed  by  the  geometrical  spreading  of  the  sound 
beam.   Thus,  for  this  simple  monostatic  situation,  we  have 


TS'  «  Sv  +  10  log  gl  +  TL^ 


+  10  log  kTT  •*   DI 
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SURFACE  REVERBERATION 


~7* 


In  the  case  of  scattering  which  arises  from  a  surface  (such  as 
the  sea  surface  or  the  bottom)  or  from  some  well-defined  layer 
such  as  the  bubble  layer  underlying  the  water-air  interface, 
the  reverberation  target  strength  is  expressed  as 


TS 


So  +  10  log  A 


where  A  is  the  surface  area  generating  reverberation  in  competition 
with  the  desired  target  and  Ss  is  the  scattering  strength  (per 
unit  area) . 

As  before,  the  problem  is  one  primarily  of  geometry ♦  involving 
the  way  in  which  source  and  receiver  major  lobes  overlap.   If  we 
again  restrict  ourselves  to  the  simple  monostatic  geometry  used 
previously  (for  which  the  major  lobe  of  the  receiver  completely 
overlaps  that  of  the  source),  then  the  surface  area  of  interest 
can  be  seen  to  be 

A  =  h   r  6  c  t 


where  9  is  the  angle  in  radians  of  the  horizontal  angular  width 
of  the  radiation  pattern  of  the  source.   This  then  gives  the 
equation 


TS*  *  S  +10  log  gl  +  10  log  r  +  10  log 
s  ? 


9 
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REPRESENTATIVE  SCATTERING  STRENGTHS 

The  values  quoted  here  are  representative  only.   Large 
deviations  can  occur  depending  on  location,  season,  biological 
activity,  and  so  forth. 

Deep  Scattering  Layers:   The  scattering  strength  Sv 
can  vary  between  about  -90  to  -60  dB  within  the 
various  sublayers. 

Water  Volume:   In  the  clear  open  ocean  waters,  Sy 
can  be  found  to  range  between  -100  to  -?0  dB. 

Surface  Scattering:   For  reasonably  low  grazing  angles, 
the  combined  scattering  from  the  ocean  surface  and 
the  population  of  bubbles  found  at  shallow  depths 
below  the  surface  gives  scattering  strengths  Ss 
ranging  around  -50  to  -30  dB  for  sea  states 
between  1  and  k. 

Bottom  Scattering:   For  grazing  angles  between  20° 

and  60°,  representative  values  of  the  scattering 
strength  Ss  range  between  -^4-0  and  -10  dB.   The 
exact  values  depend  very  strongly  on  the  particular 
bottom  and  its  composition.   In  general,  however, 
as  the  grazing  angle  goes  to  zero,  the  scattering 
becomes  negligible  and  Ss  goes  to  - oo  dB. 

REVERBERATION  LEVEL  DEPENDENCIES 

Notice  that  for  either  volume  or  surface  reverberation 
the  reverberation  target  strength  TS*  increases  with  increasing 
range,  either  as  TL„  (volume)  or  10  log  r  (surface).   This  is  to 
be  compared  with  amoient  noise  masking,  for  which  there  is  no 
range  dependence.  This  means  that,  given  sufficient  source  level, 
as  SL  M increases SONAR  performance  will  change  from  ambient- 
noise  limited  performance  to  reverberation  limited  performance. 
While  detection  for  ambient  noise  limited  performance  depends 
on  the  source  level  and  directivity  index,  notice  that  when  the 
performance  becomes  reverberation  limited  there  is  no  dependence 
on  the  source  level:   additional  power  will  not  improve 
reverberation  limited  performance.  Reverberation  can  be  reduced 
by  increasing  the  directivity  of  the  source  and/or  receiver,  or 
by  decreasing  the  pulse  length,  given  by  ct.   Increased  directivity 
has  the  disadvantage  of  reducing  scanning  rates,  and  decreased 
pulse  length  may  result  in  reduction  of  the  target  strength  TS 
of  the  desired  target. 


104 


DOPPLER  SHIFT  FOR  ACTIVE  SONAR 

In  the  case  of  active  SONAR,  any  motion  of  the  platform 
containing  the  source  will  cause  the  signal  propagated  into 
the  water  to  have  a  frequency  shift  which  depends  on  the 
bearing  into  which  the  signal  is  sent.   This  doppler-shifted 
signal  is  subjected  to  a  further  doppler-shift  when  it  is 
reflected  from  the  moving  target.   Without  going  into  the  specifics 
of  the  various  frequency  shifts,  we  shall  simply  quote  some  results 
of  interest: 

Because  of  the  two-way  propagation  of  acoustical  energy  from 
the  moving  source  to  the  moving  target  and  then  reflection  back 
to  the  moving  receiver,  the  required  bandwidth  for  the  receiver 
of  an  active  sonar  operated  under  monostatic  conditions  is 

w  in  Hz  «  1.4* (Range  rate  in  knots)* (Frequency  in  kHz) 

Notice  that  this  is  twice  the  bandwidth  required  in  the  passive 
SONAR  case. 

If  the  target  is  stationary  with  respect  to  the  surrounding 
water,  then  the  desired  echo  will  have  the  same  frequency  as  the 
reverberation  and  the  ability  of  the  receiving  system  to  recover 
signal  from  noise  will  be  minimized.   On  the  other  hand,  when  the 
target  is  moving  through  the  water,  then  the  echo  and  reverberation 
will  have  different  frequencies  and  detection  will  be  much  easier. 
This  effect  is  quite  easy  to  observe  by  ear»  The  detectibility 
of  an  "up-doppler"  or  "down-doppler"  echo  in  reverberation  is 
very  easy  even  when  the  relative  loudness  of  the  desired  echo  is 
rather  weak.   When  there  is  no  frequency  change  between  echo  and 
reverberation,  however,  the  ability  of  the  ear  to  detect  the 
echo  is  drastically  reduced.   The  strong  increase  of  the  detection 
threshold  for  slow  targets  points  to  the  use  of  parallel  filters 
of  relatively  narrow  bandwidths  to  facilitate  the  separation  of 
signal  from  reverberation  when  the  difference  in  frequency  is 
small. 
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SL     -      2TL     +     TS      ^      ML     +    (DT 
DETECTION  THRESHOLD 


Energy  Detection 

If,  as  was  the  case  for  the  passive  SONAR  equations,  the 
mechanism  for  the  detection  of  the  desired  signal  is  based  on 
detecting  an  excess  of  acoustical  energy  in  the  frequency  band 
of  interest,  then  we  have  the  same  expression  for  the  detection 
threshold, 

DT  »  5  log  ^ 

where  t  is  the  integration  time.   In  the  case  of  active  sonar, 
t  should  be  the  length  of  the  received  echo.   If  the  integration 
time  is  less  than  or  greater  than  this  pulse  duration,  it  can  be 
shown  that  the  DT  is  increased  above  this  expression,  and 
performance  degraded. 

Correlation  Detection 

An  alternative  mode  of  signal  processing  is  possible  in 
the  case  of  active  SONAR.   Since  the  amplitude  and  frequency 
properties  of  the  tone  burst  generated  by  the  source  are  known, 
it  is  possible  to  search  for  a  signal  of  these  same  properties 
in  the  received  echo.   If  the  detailed  shape  of  the  received 
echo  matches  that  of  the  sent  pulse,  then  it  can  be  shown  that 
the  detection  threshold  is  given  by 

DT  ■  10  log 


2wt 


when  the  processing  is  that  of  correlation  of  the  received  signal 
against  a  replica  of  the  generated  pulse.   There  are  many 
modifications  of  this  basic  idea,  but  all  rely  on  the  technique 
of  multiplying  the  received  signal  by  a  time-delayed  model  of 
the  sent  pulse  and  integrating  the  resultant  product  over  the 
pulse  duration.   This  is  done  for  monotonically  increasing 
delay  times,  and  the  resulting  function  of  delay  time  studied 
for  the  presence  of  a  sharp  peak  representing  the  matching  of 
the  delayed  replica  with  the  same  signal  appearing  in  the  received 
signal  and  noise.   In  practice,  the  finite  size  of  the  target, 
multipath  interference,  and  the  fluctuations  introduced  by  the 
inhomogeneities  of  the  transmission  of  the  signal  through  the 
water  alter  the  received  echo  so  that  the  correlation  process 
is  not  as  good  as  would  be  desired.   As  a  result,  the  DT  for 
a  real  system  operating  in  the  real  world  is  more  than  the 
theoretical  expression  would  predict. 
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SUMMARY  OF  THE  ACTIVE  SONAR  EQUATIONS 


Ambient  noise  limited  conditions 


SL  -  2TL  +  TS  ^  HBL  -  AG  ♦  10  log  w  +  DT 


Reverberation  limited  conditions 


TS  -  TS  ■  >  DT 


<? 
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EXAMPLE 

An  active  SONAR  operates  at  1  kHz  with  a  source  level  of  220  dB  re  1  fPa, 
a  directivity  index  of  20  dB,  a  horizontal  beam  width  of  10°,  and  a  pulse 
length  of  0.1  sec.  Correlation  detection  is  used  and  it  is  desired  to  have 
P(D)  -  e.50  with  P(FA)  -  10-^.  The  target  strength  of  the  submarine  is 
expected  to  be  30  dB  and  its  speed  may  be  up  to  20  kts.  Both  the  SONAR 
and  the  submarine  are  in  a  mixed  layer  of  depth  100  m  with  the  source  at 
36  m.  The  sea  state  is  3>  the  scattering  strength  for  surface  scatterers 
is  -30  dB,  and  volume  scattering  is  negligible.  For  simplicity,  assume  that 
conditions  are  such  that  the  detection  threshold  is  the  same  for  both  noise- 
limited  and  reverberation-limited  performance. 
Solution: 

w  ■  l,k   (frequency  in  kHz)  (range  rate  in  kts)  ■  l,k  x  1  x  20  ■  28  Hz 

d  ■  15      same  calculation  as  in  example  1  (passive) 
UT  -  10  log  (d/2wt)  -  10  log  (15/2x26x0.1)  -  10  log  2.68  -  5  dB 
NL  -  NSL  -  IE  +  16  log  w  -  62  -20  +  LU  -  56  dB 
r+  ■  1310  m     and    r  ■  8u00  m     same  calculation  as  Example  3  (passive 

Assume  that  the  detection  is  noise  limited 

2TL  ^  SL  +  TS  -  NL  -  DT  -  220+30-56-5  .  189 
TL  ^  95 


-4 


TL  -  31  +  10  log  r  +2.6  x  10 

10  log  r  +  2.6  x  10"4  r  =  95-31  =  64 


same  calculation  as  Example  3 


r 

10 

log  r 

!04 

40 

105 

50 

4 
5   x  10 

47 

6  x   10 

48 

-4 
2.6  x  10   r 


r  =  6  x  10  m 


26 
26 
13 
16 

if  noise  limited 


TL 


43 
76 
60 
64 
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EXAMPLE  (CONTINUED) 

But,  vhat  is  the  RL  at  this  range? 

TS»  »  Ss  +  10  log  (ct/2)  +  10  log  r  +  10  log  6 

-30    +     10  log  (1^00  x  0.1/2)     +48      +10  log  0.17U 
-30+19+48-8     -     29    dB 
RL       -       SL     -     2TL     +     TS« 

220  -  189  +    29     .    60    dfi 
So 

RL  S  NL 
Therefore  the  detection  will  be  reverberation  limited 
TS  -  TS  •    >,    DT 

TS'     -     -30  +  19  +     10  log  r    -  8      -      -19    +     10  log  r 
TS     -    TS'    >/    m 
30+19     -     16  log  r      >,     5 
10  log  r    ^    hk 

r     -       2.5  x  10^  m 
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APPENDIX  A:   LOGS 
All  you  have  ever  wanted  to  know  about  the  decibel  but  were  afraid  to  ask, 
The  language  in  which  the  SONAR  equation  is  written  is  the  language  of  the 
decibel.    This  self-paced  exercise  is  intended  to  make  you  familiar  with  the 
grammar  of  this  language. 

Instructions.   For  each  numbered  section,  first  attempt  to  do  the  problems. 
(The  answers  are  on  the  last  page.)   If  you  have  no  troubles,  skip  over  to  the 
next  section.   On  the  other  hand,  if  you  have  difficulties^  read  through  the 
discussion  and  then  try  the  problems  again.   Go  to  the  next  section  only  when 
you  can  do  the  problems  without  difficulty.   By  the  time  you  finish  all  the 
sections,  you  should  be  speaking  "decibel"  like  a  native. 

1.   EXPONENTIAL  NOTATION 
Problems  A 

1.  Write  each  of  the  following  numbers  in  exponential  forms 

a)  2,365,000 

b)  872 

c)  3340 

d)  42.1 

e)  600 

f)  1.0 

2.  Write  each  of  the  following  numbers  in  conventional  form 

a)  2.689  x  102 

b)  7.30  x  102 

c)  6.29  x  104 

d)  8.9  x  10 

e)  3.261  x  102 

f)  1.00  x  10° 
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3.  Write  each  of  the  following  numbers  in  exponential  form 

a)  0.00252 

b)  0.012 

c)  0.00002 

d)  0.01002 

e)  0.00328 

4.  Write  each  of  the  following  numbers  in  conventional  form 

a)  1.03  x  10~6 

b)  3.2  x  10"2 

c)  6.89  x  10~4 

d)  1.003  x  10"1 

e)  1.00  x  10"1 

A  succinct  way  of  expressing  very  large  and  very  small  numbers  is  by 
exponential  (or  scientific)  notation.   For  example,  one-hundred  and  twenty-six 
million  can  be  written  either  as 

126,000,000 
or  as 

1.26  x  108  , 
where  the  exponent  8  denotes  the  number  of  places  that  the  decimal  must  be 
moved  to  the  right  to  express  the  number  in  its  conventional  form. 

There  is  no  unique  way  of  expressing  a  number  in  exponential  form,  since 
1.26  x  108  =  12.6  x  107  =  0.126  x  109  =  etc. 
Which  form  you  use  depends  on  convenience. 

Numbers  smaller  than  unity  can  be  handled  in  a  similar  manner.   For 
example,  one-hundred  and  twenty  millionths  can  be  written  either  as 

0.000120 


or  as 


1.2  x  10~4  , 
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where  the  negative  exponent  denotes  that  the  decimal  place  is  to  be  moved  to 
the  left. 

A  special  case  that  must  be  watched  is 

1  -  1  x  10°  -  10° 
which  is  consistent  with  the  above  conventions. 
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2.   ALGEBRA  OF  EXPONENTS 
Problems  B 

1.  Calculate  the  following  and  express  the  answer  in  exponential  form 

a)  6.15  x  103  +  2.34  x  103 

b)  3.2  x  102  -  1.46  x  103 

c)  9.1  x  104  +  2.4  x  105 

d)  3.58  x  10~2  +  1.26  x  lo"2 

e)  6  x  10~6  +  3.281  x  10~4 

f)  3.5  x  104  -  2  x  101 

2.  Calculate  the  following  and  express  the  answers  in  exponential  form 

a)  2.50  x  106  x  2.00  x  10 

b)  4.40  x  104  /  2.20  x  105 

c)  3.26  x  103  x  2.00  x  10~2 

d)  5.34  x  10~2  /  2.67  x  105 

e)  8.68  x  10_1  /  4.34  x  lo"2 

3.  Calculate  the  following  and  express  the  answer  in  exponential  form 

a)  (1.1  x  104)2 

b)  (2.0  x  10~3)4 

c)  (1.44  x  102)1/2 

d)  (1.44  x  10_2)1/2 

e)  (8.0  x  106)1/3 

4  0.5 

f)  (2.56  x  10  ) 

x   /o    ,  -6,0.3333 

g)  (8  x  10   ) 

To  add  and  subtract  numbers  expressed  in  exponential  form  it  is  necessary 

4 

to  express  both  numbers  in  the  same  exponent.   For  example,  adding  3.56  x  10 

to  2.1  x  10  is  accomplished  by 

3.56  x  104 
0.21  x  104 

3.77  x  104 
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-1  3 

When  asked  to  subtract  8.67  x  10  "  from  2.23  x  10  ,  note  that  the  answer  is 

2.23  x  10  ,  for  the  term  to  be  subtracted  affects  the  larger  term  beyond  the 

figures  retained. 

Exponential  notation  really  shows  its  worth  in  calculations  involving 

multiplication  and  division.   We  see  that 

120  x  2000  -  240,000 

or  1.2  x  102  x  2  x  103  =  2.4  x  105 


and 


or 


36000  /  200  =  180 


3.64  x  104  /  2  x  102  =  1.8  x  102. 
The  rules  are  simple: 

In  multiplication,  exponents  add. 
In  division,  exponents  subtract. 

A  special  case  that  is  often  useful  is 

2      0      2-2 
1  /  10   =  10U  /  10^  =  10  *, 

which  allows  exponents  to  be  moved  between  numerator  and  denominator. 

Power  and  roots  are  trivial  in  exponential  notation: 

(1.2  x  103)2  =  1.2  x  103  x  1.2  x  103  =  1.44  x  106, 


or  for  any  power  n 


and, 


or  for  any  power  of  n 


(10V  =  10™ 


(1.44xl0V/2=  (1.44)1/2x  (106)1/2  =  L2  X  103' 


(ioV/n  =  iom/n. 


This  latter  form  introduces  the  possibility  of  fractional  exponents;  these  will 
be  discussed  in  Section  4. 
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3.   LOGS  (IN  BASE  TEN) 
Problems  C 

1.   Find  x 

3 

a)  log  10  =  x 

b)  log  10"6  =  x 

c)  log  10  =  4 

d)  log  10X  =  -3 

e)  log  10  =  0 

Since  we  say  that  "ten  to  the  fourth  power  equals  ten  thousand"  we  can 
ask  "to  what  power  must  ten  be  raised  to  get  ten  thousand?"  The  answer  is 
obviously  "four."  You  are  now  an  expert  on  logarithms,  since  the  "log"  of 

10,000  is  4.   This  is  written 

4 
log  10,000  =  4  or  log  10  =  4 

and  this  is  all  there  is  to  logs. 

Find  the  log  of  x  by  asking  yourself  to  what  power  must  I  raise  10  to 
get  x. 

It  is  equally  appropriate  to  ask  "the  log  of  what  is  4,"  that  is  what 
is  x  so  that  log  10  =4. 
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4.  DECIMAL  EXPONENTS 


PROBLEMS 

D 

1. 

Find 

a) 

log 

4 

b) 

log 

40 

c) 

log 

0.4 
.2 

d)  log  4 
e) "  log  1/4 
2.   Find  the  number  whose  log  is 

a)  0.70 

b)  1.70 

c)  -1.30 

d)  -1.70 

e)  -2.70 

So  far  we  have  assumed  integer  exponents,  but  everything  we  have 
said  also  is  true  for  decimal  exponents.   For  example, 

10°-3°=2.0 
and 

log  2  =  0. 30 
The  table  below  will  allow  you  to  calculate,  with  a  fair  degree  of  accuracy, 
all  logs. 

1.0    1.2    1.5    1.7    2.0   .3.0    4.0    5.0    6.0    7.0    8.0   9.0 


10.0 


log  x 


0.00   0.08   0.18   0.23   0.30   0.48  0.60   0.70   0.78   0.85   0.90  0.95   1.0< 
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If  you  ever  need  more  accuracy  you  will  have  to  use  slide  rule,  calculator, 
or  tables. 

To  find  the  log  of  a  number  greater  than  10  proceed  as  follows: 
log  7000  =  log  7  x  103  =  log  7  +  log  103 
=  0.85  +3    =  3.85 


or 

log  0.003=  log  3  x  10~3=  log  3  +  log  10~3 
=  0.48  -  3    =  -2.52 
There  is  one  property  of  logs  that,  if  used  correctly,  will  greatly 

simplify  many  calculations.   This  will  be  illustrated  by  an  example: 

2 
log  36  =  log  6   =  log  (6x6)  =  log  6  +  log  6=2  log  6=  1.56 

(  Calculate  log  36  and  compare  to  2  log  6 . ) 

In  general 

,    n 

log  x  =  n  xog  x. 

It  is   also   true   that 

log  xy  =   log  x  +   log  y. 

log  36  =  log  (4x9)  =  log  4+2  log  3  =  0.60  +  0.96  =  1.56 
But,  remember  that 

log  (x  +  y)  is  not  equal  to  log  x  +  log  y  ! ! ! 
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5.  THE  DECIBEL 

PROBLEMS  E 

1.  Express  these  signal-to-noise  ratios  in  dB 

a)  2 

b)  1/2 

c)  10 

d)  1/10 

e)  100 

2.  The  voltage  of  a  signal  is  1  volt*,  express  the  signal  level  with 

reference  to  the  following  voltages 
.-12 


a) 

10 

b) 

io-6 

c) 

1 

d) 

io6 

e) 

io12 

By  convention,  when  the  log  of  a  pressure  or  voltage  is  multiplied 
by  20,  the  resulting  number  is  said  to  be  measured  in  "decibels",  donated 
by  dB. 

20  log  2  =  6dB 
20  log  10  =  20dB 
Decibels  are  used  in  two  different  ways: 

a)  to  express  a  ratio  such  as  the  signal-to-noise  ratio 

signal- to- noise  ratio  =  20  log  (signal  voltage/noise  voltage) 

b)  to  express  a  level  such  as  signal  level 

signal  level  =  20  log  (signal  voltage/reference  voltage) 
In  the  latter  case  it  is  absolutely  essential  to  specify  the  reference, 
e.g.,  "the  signal  level  is  120dB  re  10   volts". 
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ANSWERS 
Problems  A:    l.a)      2.365  x  106,   b)    8.72   x  102,    c)    3.34   x  103,    d)    4.21   x  101, 

e)  6   x  102,    f)    1.0   x  10° 

2. a)  268.9,    b)    730,    c)    62,900,    d)    8,900,000,    e)    326.1 

f)  1.00 

3. a)  2.52   x  lO-3,   b)    1.2   x  10"2,    c) 2  x  10"5,    d)    1.002   x  10~2 

e)  3.28  x  10~    . 

4. a)  0.00000103,    b)    0.032,    c)    0.000689,    d)    0.1003,    e)    0.100 


Problems  B:    l.a)    8.49   x  103,   b)    -1.14   x  103,    c)    3.3  x  105,    d)    4.84  x  10~2 
e)    3.34  x  10"4,    f)    3.5   x  104. 

2. a)    5.00  x  1010,   b)    2.00   lO-1,    c)    6.52   x  101,    d)    2.00   x  lO-7, 
e)    2.00  x  10. 

3. a)    1.2   x  108,   b)    1.6   x  10-11,    c)    1.2  x  101,    d)    1.2   x  lO-1, 
e)    2.0   x  102,    f)    1.6  x  102,    g)    2.    x  10~2. 


Problems  C:    l.a)  3,   b)      -6    ,c)    4,    d)    -3,    e)    0 

Problems   D:    l.a)  0.60,    b)    1.60,    c)    -0.40,    d)    1.20,    e)    -0.60. 

2. a)  5.0,   b)    50,    c)    0.050,    d)    0.020,    e)0.0020 

Problems   E:    l.a)  6dB,    b)    -6dB,    c)    20dB,    d)    -20dB,    e)    40dB. 

2. a)  240dB,   b)    120dB,    c)    OdB,    d)    -120dB,    e)    -240dB 
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